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ABSTRACT 


The  proposed  reliability-based  static  strength  criteria  system  described  in  AFFDL- 
TR-67-107,  Volumes  l-lll,  was  reviewed  to  determine  the  data  requirements 
and  availabi.ity,  the  implications  of  such  an  approach  on  the  structural  design 
process,  methods  by  which  implementation  can  be  achieved  without  discontinuity, 
and  necessary  changes  to  specification  and  handbooks.  Volume  I  describes  the 
studies  made  using  data  for  the  C— 141  cargo  transport.  Volume  II  describes  the 
findings  and  includes  five  appendices.  The  principal  conclusions  are  that  insufficient 
data  exists  for  the  imminent  implementation,  but  that  studies  of  the  relative  reliability 
of  different  configurations  and  components  or  of  different  conditions  at  the  same 
location  would  provide  a  short  term  means  of  using  the  system  to  gain  familiarity 
and  confide  ce. 
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SECTION  I 


INTRODUCTION 


Many  attempts  have  been  made  to  achieve  the  realization  of  techniques  for  applying 
reliability  methods  to  the  definition  of  structural  strength.  The  most  comprehensive  of 
these  was  prepared  by  Innes  Bouton  and  others  and  is  described  in  AFFDL-TR -67-107. 
The  three  volumes  of  that  report  of  that  report  discussed  previous  methods  and  derived 
proposed  methods  covering  both  time-independent  (static)  and  time-dependent  (fatigue) 
strength.  The  full  range  of  interactions  with  non -structural,  operational,  executive, 
and  contractual  areas  was  discussed. 

The  study  described  in  the  present  report  was  aimed  reviewing  the  proposed  method 
for  applying  probabilistic  techniques  to  the  assessment  of  static  strength  reliability. 

This  review  was  to  identify  the  data  requirements  of  the  proposed  method,  the 
necessary  changes  to  specifications  and  design  handbooks,  the  interfaces  with  non- 
structural  design  areas  and  the  steps  to  be  taken  during  implementation  of  the  method. 


(Repeated  from  Volume  I.) 
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SECTION  II 
SUMMARY 


A  cieor  understanding  of  the  various  operations  incorporated  into  the  proposed  static 
strength  reliability  analysis  of  AFFDl-TR-67-107  is  necessary  to  its  successful  imple¬ 
mentation  .  Section  III  provides  a  simple  worked  example  which  illustrates  each  step 
in  turn  using,  first,  dummay  date  and  then  realistic  data.  The  categories  of  required 
data  are  defined . 

Sections  IV  through  IX  discuss  each  category  in  turn,  by  means  of  studies  of  data 
pertinent  to  the  C-14IA  cargo  transport  aircraft.  Section  X  then  summarises  the 
findings  in  the  form  of  a  trial  application  of  the  method  to  the  wing  of  the  C-I4IA. 

Sections  XI  and  XII  discuss,  respectively,  the  updating  of  the  data  to  reflect  the  state 
of  knowledge  at  each  stage  during  the  design  and  operational  life  of  a  vehicle,  and 
the  form  in  which  the  required  doto  might  be  standardised. 

Specific  steps  reouired  to  achieve  the  short-term  and  long-term  implementation  of  the 
method  ore  described  in  Section  XIII,  and  the  necessary  changes  to  existing  Mll-A 
specifications  and  AFSC  Design  Handbooks  ore  summarised  in  Section  XIV.  Section 
XV  contains  the  conclusions  and  recommendations  resulting  from  the  study. 

Five  appendices  follow  the  main  text.  Appendix  i  outlines o  technique  for  the  use  of 
bi -modal  (double -fondly)  statistical  distributions;  the  Gumbo  I  distribution  of  extremes 
is  employed  as  ©n  example,  but  the  method  is  valid  for  o  range  of  statistical  distribu¬ 
tions.  Appendix  II  contains  the  basic  equations  of  the  computer  program  used  in  the 
study;  this  uses  double -family  Gumbo  I  distributions,  a  constant  calculation  interval, 
end  employs  So  yes'  theorem  to  incorporate  the  effects  of  test  results,  but  is  otherwise 
similar  to  the  origincl  program;  many  of  the  informed  iote  results  are,  however,  printed. 
Appendix  III  describes  the  program,  its  input  requirements  and  operation. 

Appendix  IV  contains  sample  runs  mode  with  the  program,  and  Appendix  V  shows  tlx; 
analysis  of  load  end  strength  data  using  double-fcmily  representations. 

(Repeated  from  Volume  !.) 
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SFCT10N  XI 

UPDATING  OF  DATA 


ll.l  Introduction 

Reference  l  stresses  the  continuous  nature  of  the  process  of  establish¬ 
ing  the  structural  reliability.  The  specific  items  to  be  updated  ore 
described  in  this  Section,  *ith  the  practical  means  of  doing  so  described. 


1 1 .2  Data  Items 

The  three  fundamentoJ  data  categories  ore: 
o  load  spectrum 

O  error  function 

o  strength  distribution 

and  each  will  change  periodically  during  the  total  lifetime  of  a  specific 
aircraft.  The  particular  points  at  which  data  revisions  ore  most  likely, 
and  which  permit  progressively  updated  reliability  estimation,  are: 
a  Initial  Design  Stage 

o  Detail  Design  Stage 

a  After  Detail  Design,  but  before  Static  Testing 

o  After  Static  Testing,  but  before  Design  Revision 

Final  Design,  but  before  Operation 
o  During  Operation 

bach  is  discussed  separately  below . 


1 1 .3  Revision  $1 


initio!  Design 

During  initial  design,  the  load  spectra  must  be  based  on  assumed 
utilisation,  assumed  aerodynamic  and  inertia  distributions  over 
the  airframe  and  assumed  probabilities  of  occurrence  of  different 
conditions.  An  error  function  con  offeitrorily  be  selected  from  one 
of  tne  "standard  set",  or  can  be  based  on  post  test  experience 
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within  the  particular  company.  Strength  distribution  data  will 
be  selected  from  the  chosen  material  data,  with,  advisedly, 
allowances  for  the  effects  of  fabrication  and  assembly  which 
reflect  any  unconventional  feotures .  Predictions  can  be  mode 
of  the  reliability,  assuming  values  for  the  various  parameters. 

b.  Detail  Design 

By  the  time  thot  the  detail  design  stage  is  reached,  some  additional 
information  will  generally  be  available  .  Revised  load  spectra  will 
have  replaced  the  preliminary  data;  some  component  test  data  will 
usually  have  been  accumulated,  particularly  for  any  novel  design 
features,  and  will  permit  a  revised  error  function  to  be  selected. 

If  new  construction  methods  are  proposed  (fasteners,  soy',  then 
sufficient  rest  data  will  perhaps  be  available  to  indicate  the 
variability  of  the  process  and  so  to  permit  revision  of  the  strength 
distribution,  A  second  set  of  reliability  estimates  is  possible . 

c.  Before  Static  Testing 

At  the  end  of  the  detail  design  stage,  but  before  static  testing, 
a  third  vet  of  reliability  estimates  con  he  calculated*  This 
w»|f  reflect  any  additional  data  gathered  up  to  this  time,  par¬ 
ticularly  In  the  strength  distribution  orec.  The  reliability 
prediction*  will  re.noin  hosed  on  assumed  test  results . 

d.  After  Static  Testing 

The  static  test  results  will  Have  one  of  two  effects.  Either 
the  dbsign  goal  will  hove  been  met,  th  confirming  the  pre¬ 
dictions,  or  it  will  not  hove  been  met .  |n  the  loiter  event, 
two  courses  of  oct'on  ore  possible:  redesign  will  be  performed 
in  the  foiled  regions,  representing  a  further  iteration,  or  the 
design  and  operating  conditions  will  be  revised  ta  correspond  to 
achievement  of  a  lower  loading  of  the  original  reliability,  or  a 
lower  reliability  a*  the  original  Ic-od. 
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e.  Final  Design 

After  any  redesign  or  re-analysis  has  been  completed,  but  before 
the  aircraft  enters  service,  a  further  reliability  assessment  can  be 
made.  This  will  still  be  based  on  assumed  utilization  ond  assumed 
load  distribution  data,  but  will  reflect  all  strength  data  accumulated 
up  to  this  time . 

f.  During  Operation 

Operational  data  will  be  appropriate  to  two  distinct  types  of  revision 
of  the  reliability  estimate.  The  first  is  the  obvious  one  of  permitting 
realistic  load  spectra  to  be  formulated,  ond  the  second  is  a  very 
important  one  which  is  usually  overlooked.  Each  flight  experience  of 
a  particular  load  is  on  additional  test  to  that  load  level.  Now,  it 
has  been  shown  that  the  influence  of  testing  to  low  load  Itwels  is 
insignificant ,  but  each  and  every  aircraft  that  experiences  a  high 
load  level  provides  o  further  data  point  which  adds  to  the 
knowledge  required  to  nrodiet  a  better  reliability. 

Periodic  updating  can  be  performed  os  data  is  accumulated; 
this  should  not  be  too  frequent,  for  economic  reasons,  and 
determination  of  the  appropriate  times  will  depend  cm  individual 
circumstances , 

U  .4  Operational  Pott?  Recording 

Q.  One  of  the  greatest  potential  areas  for  acquiring  new  and  better 
structural  design  data  lies  within  the  Air  Force’s  Aircraft  Struc¬ 
tural  Integriry  Program,  A$|P,  As  a  port  of  the  AS  IP,  each  aircraft 
system  must  have  an  Individual  Aircraft  Service  Life  Monitoring 
Program,  IAS  IMP,  and  as  a  port  of  the  IASIMP  for  the  more  critical 
systems,  a  number  of  aircraft  in  each  fleet  is  to  be  equipped  with 
Multi- Channel  fiecorden  (M-CRK 
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The  need  for  multi-channel  recorders  has  been  a  recognized  part  of 
Air  Force  planning  for  at  least  ten  years .  Some  recorders  have  even 
been  developed  end  used  with  varying  degrees  of  success,  but  with 
limited  applicability.  Starting  in  1968,  the  Air  Force  laid  plans 
for  a  new  and  more  universal  recording  system.  The  AFLC,  through 
its  several  AMA's,  gathered  data  on  the  type  of  information  needed 
to  effectively  carry  out  the  lASLMP's  on  a  wide  variety  of  aircraft. 
These  data  were  synthesized  by  ASD,  along  with  other  known  and 
projected  requirements,  to  prepare  a  set  of  recorder  specifications. 

In  June  1970  ASD,  under  the  auspices  of  AFLC,  let  a  contract  to 
develop  a  new  24  channel  digital  recording  system  and  a  ground 
playback  unit .  That  system  is  still  under  development  at  this  writing  . 
A  unique  feature  of  the  new  system  is  that  a  single  basic  recorder 
unitwill  be  suitable  for  all  types  of  aircraft.  To  accommodate 
peculiar  requirements  of  different  types  of  aircraft,  the  system  in¬ 
cludes  development  of  four  different  converter/multiplexer  units, 
each  of  which  is  compatible  with  the  one  recorder  module.  Current 
plans  call  for  the  initial  production  of  about  140  recorder  systems, 
with  a  contingency  buy  ot  approximately  140  additional  S) stems. 

A  portion  of  most  "first-line"  aircraft  fleets  (ranging  from  about 
5  to  20%)  are  tentatively  scheduled  to  receive  the  recorders, 
wirh  t«  st  installation  starting  in  late  1972. 

L.  One  of  the  major  objectives  of  the  multi-channel  recorder 
program  is  to  provide  a  better  tool  by  which  to  accomplish 
structural  fatigue  tracking.  In  fact  the  entire  program  to 
date  has  been  oriented  toward  -  and  largely  justified  by  - 
the  structural  fatigue  problems.  However,  because  of  the 
htgh  commonality  between  the  data  needed  for  fatigue  design 
or  tracking  and  the  data  required  to  develop  new  statistically 
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TABLE  XIX 


PROPOSED  LIST  OF  MULTI-CHANNEL  RECORDER 


PARAMETERS  FOR  THE  C-141. 


ITEM 

NAME 

T 

Clock  Time 

HP 

Pressure  Altitude 

V 

e 

Equivalent  Airspeed 

Nz 

Normal  Acceleration  at  C.G.,  g's 

Ny 

Lateral  Acceleration  at  C.G.,  g's 

6 

Pitch  Rate 

V 

Yaw  Rate 

Elevator  Position 

Rudder  Position 

* 

Flap  Position 

Vg 

Ground  Speed 

ft, 

Nose  Gear  Steering  Angle 

*1 

Strain  at  Location  1 

°2 

Strain  at  Location  2 

«3 

Strain  at  Location  3 

®4 

Strain  at  Location  4 

°5 

Strain  at  Location  5 

AP 

Cabin  Pressure  Differential 

Wf 

Total  Weight  of  Fuel 

s.s. 

Squat  Switch  Make-or- Break  Signal 

DD1 

Date 

DD2 

Serial  Number 

DD3 

Base  of  Assignment 

DD4 

Initial  Cargo  Weight  or  Cargo  Update 

DD5 

Total  Initial  Fuel  Weight 

bosed strength  design  criteria,  this  letter  area  will  inevitably 
benefit.  A  second  stated  objective  of  the  Multi-channel  recorder 
program  is  to  accumulate  data  for  use  in  the  structural  design  of 
future  aircraft  systems ,  The  opportunity  offered  by  this  objective 
is  obvious.  After  the  recorder  installations  are  made  and  records 
accumulated  for  a  time ,  certain  of  the  data  on  each  aircraft  system 
will  "mature"  to  the  point  that  further  recording  produces  no  new 
intelligence.  When  this  happens  the  recorder  program  can  and 
should  be  redirected  toward  goals  that  are  allied  more  specifically 
to  the  statistical  strength  criteria  .  A  hypotherica!  example  of  such 
a  switch  is  given  in  the  second  following  paragraph. 

c.  Since  the  C-141  has  been  used  for  illustration  in  other  sections  of 
this  report,  the  same  theme  is  followed  here  .  Table  XIX  shows  the 
list  of  parameters  that  have  been  selected  by  WRAMA  for  the  M-CR 
program  on  the  C- 141  .  Each  of  these  parameters  -  either  singly  or 
in  combination  with  others  -  is  believed  to  produce  something  of 
value  either  directly  in  a  contemporary  fatigue  tracking  process, 
or  in  the  derivation  of  new  fatigue  criteria,  or  both.  Several  of 
these  parameters  should  be  useful  also  in  deriving  new  statistical 
strength  design  criteria.  For  example,  normal  acceleration  at  the 
center  of  gravity,  N^,,  is  a  valuable  parameter  in  its  own  right 
since  it  is  a  direct  measure  of  the  gross  symmetric  response.  When 
experience  is  properly  sorted  by  weight,  speed  and  altitude, 
subsequent  peak  counting  of  the  data  provides  the  type  of  statistical 
distributions  required  in  the  determination  of  strength  levels,  A 
further  refinement  is  attainable  by  a  joint  analysis  of  with 
elevatui  deflection.  Besides  allowing  a  separation  of  symmetric 
gust  and  maneuver  responses,  such  an  analysis  should  also  afford 
the  collection  of  good  statistical  samples  of  abrupt  pitch  maneuvers  - 
an  important  design  area  about  which  little  is  known  for  cargo  aircraft. 
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In  a  similar  manner  Ny  can  be  sorted  and  peak  counted  for  use 
in  lateral  load  predictions,  and  it  can  be  jointly  analyzed 
with  rudder  deflection  to  fill  the  gap  in  knowledge  about  abrupt 
rudder  kick  maneuvers . 

d.  One  other  illustrative  example  is  worthy  of  noting  here  concerning 
possible  future  direction  of  the  M-CR  program  .  If  the  C- 141 
program  is  successfully  implemented  and  prosecuted,  certain  of  the 
parameters  will  attain  a  statistical  stability  after  which  they  need 
not  be  recorded  full  time.  The  resulting  surplus  of  recorder  capacity 
may  be  used  effectively  to  fill  knowledge  gaps  such  as  that  concerning 
the  phasing  of  PSD  loads.  In  gust  analysis  current  PSD  methods  allow, 
for  example,  a  fairly  precise,  but  separate,  determination  of  shear, 
bending  and  torsion  at  a  given  structural  location;  but  the  phasing 

of  these  three  vectors  in  a  deterministic  strength  analysis  is  largely 
a  guessing  game  .  The  addition  of  strain  gage  clusters  or  rosettes 
at  selected  stations  could  provide  real  It'-'  samples  of  the  amplitude 
and  frequency  relationships  among  two  or  more  load  vectors.  Data 
such  as  these  will  be  essential  in  future  designs  in  order  to  express 
applied  loads  and  structural  strength  in  a  common  set  of  terms. 

e.  In  summary,  the  multi  -channel  recorder  program (s)  is  viewed  as 
havinr,  great  potential  in  gathering  data  for  use  in  applying  the 
new  structural  reliability  concepts.  In  particular,  statistical 
information  on  loads  and  the  loads  environment  will  be  developed 
in  both  the  volume  and  detail  necessary  for  a  rigorous  statistical 
load  analysis . 
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SECTION  XII 
STANDARDIZED  DATA 


12.1  Data  Categories 

a.  The  proposed  probabilistic  system  of  criteria  is  intended  to 
provide  desired  structural  reliability  levels  for  normal  operation 
and  for  a  reasonable  degree  of  overload.  The  basic  concept  is 
that  the  statistical  variations  of  load  and  strength  are  jointly  in 
corporated  into  the  risk  assessment.  A  convenient  approach  would 
be  provided  by  charts  relating  the  required  design  and  test  factors 
to  the  reliability  levels  in  terms  of  parameters  describing  the  load 
strength  distributions,  the  error  function  and  the  number  and  type 
of  tests. 

b.  The  study  described  in  this  report  reveals  that  the  "demonstrated" 
reliability  is  indeed  a  function  of  all  of  these  parameters.  Further¬ 
more,  simultaneous  consideration  of  both  normal  (limit)  and  overload 
(omega)  conditions  will  seldom  be  possible  because  the  permissible 
load  level  (streng  will  generally  be  different.  Each  structural 
location  will  require  separate  analysis,  due  to  the  fact  that  both 
the  load  and  stre.  gth  distribution  will  differ  from  point  to  point. 

12.2  Load  Data 

a.  Theoretically,  it  would  appear  possible  to  develop  o  single  load 
spectrum  for  each  iocotion,  which  would  contain  the  total  load 
Occurrence  properties  for  on  aircraft  lifetime.  The  statistics  required 
♦o  achieve  this  goal  are  not  available,  even  on  aircraft  which  have 
accumulated  extensive  operational  experience.  For  example,  in¬ 
formation  is  required  on  the  probabilities  of 

1)  weight  and  weight  distribution 

2)  speed  and  height 

3)  type  of  load  condition  (gust,  pull-up,  rudder 
kick,  etc .) 

4)  level  of  loading  (in  terms  of  a  basic  parameter) 

5)  time  history  of  loading  (to  describe  the  local 
loading) 

6)  associated  load  systems  (pressure,  thermal  gradient, 

etc .) 

10 


and  these  probabilities  are  clearly  not  independent,  so  that 
the  resultant  probability  of  each  combination  is  needed. 

b.  Some  degree  of  standardization  may  be  feasible,  even  if  it 
is  more  arbitrary  than  statistical  in  origin.  For  example, 
gust  velocity  descriptions  are  already  employed  in  fatigue 
analysis  and  would  be  directly  usable.  Normal  load  factor 
spectra  exist  in  a  suitable  form  in  existing  criteria  (reference  2) 
and  these  distributions  can  be  regarded  as  standard  for  the 
appropriate  category  of  aircraft  and  the  appropriate  type  of 
mission. 

c.  Many  of  the  remaining  areas  require  extensive  data  collection 
and  analysis.  This  is  particularly  true  of  the  asymmetric  flight 
conditions  whit.,  are  of  increasing  significance  cr  sweepback 
increases  and  aspect  ratio  reduces. 

d.  For  the  initial  use  of  the  proposed  system,  one  possible  means 
of  filling  the  void  would  be  for  an  assumed  set  of  data  to  be 
derived  from  what  data  can  be  assembled.  Such  synthetic 
"statistics'*  must  be  regarded  os  artificial,  but  would  at  least 
permit  comparison  of  different  aircraft,  different  locations  on 
the  same  aircraft  or  different  structural  designs  of  the  same  loca¬ 
tion. 

e.  The  necessity  for  an  adequate  probabilistic  prediction  of  the 
utilization  of  the  aircraft  becomes  os  great  os  in  fatigue  analysis. 
However,  in  addition  to  the  overage  or  typical  conditions  for 
each  segment  of  the  mission  profile,  it  will  be  necessory  to 
derive  (or  ossume)  the  shape  and  dispersion  about  this  mean. 

Without  this  detailed  iovel  of  data,  no  realistic  estimate  of 

the  risk  of  failure  can  be  made;  the  alternative  is  to  ignore 
the  probability  distribution  of  the  loads,  to  assume  a  known 
certain  load  and  to  base  the  reliability  estimate  solely  on  the 
variation  in  strength.  While  this  may  be  conservative,  it  negates 
most  of  the  advantages  implicit  in  the  proposed  method  of  reference  1 


12.3  Strength  Doto 

a.  Extensive  materiol  strength  data  already  exists  as  the  necessary 
means  of  establishing  the  published  design  allowables.  Although 
the  allowables  themselves  represent  only  one  discrete  point  in 
the  distribution,  the  required  data  should  be  accessible;  the  form 
of  data  required  consists  of  the  mean  and  standard  deviation  and 
the  shape  of  the  distribution  to  be  used.  For  reliability  analyses, 
it  should  be  remembered  that  the  lower  tail  of  the  distribution  is 
most  important  in  the  assessment  of  the  risk  of  first  failure  (mean 
time  before  failure  estimates  are  not  appropriate),  and  distributions 
must  be  chosen  with  this  in  mind.  Double-family  distributions  may 
be  appropriate  (see  Appendix  I),  ond  if  these  are  employed,  the 
necessary  material  strength  data  will  generally  contain  five  parameters. 

b.  Most  of  the  data  described  above  relates  to  the  basic  properties  of 
the  material  as  delivered  to  the  aircraft  manufacturer.  The  structural 
strength  of  the  final  product  will  reflect  variations  imposed  by  all  of 
the  operations  inherent  in  fabrication  and  assembly  (the  time -dependent 
effect  of  service  wear  and  tear  is  not  considered  in  the  contoxt  of  the 
present  study,  but  may  need  to  be  examined). 

Data  on  the  strength  of  various  detail  configurations,  such  as  lugs, 
fittings,  joints,  etc.  exists  in  o  random  manner,  but  usually  in  in¬ 
sufficient  quantity  to  provide  adequate  statistical  distribution  dato. 

The  acquisition  of  such  information  is  of  paramount  importance  to 
the  success  of  the  proposed  method. 

c.  Appendix  V  gives  examples  of  the  analysis  of  typical  samples  of  dato 
of  materiol  strength  ana’  joints.  These  reveal  that  the  convention^ 
assumption  of  normal  distributions  moy  not  be  desirable,  end  that 
better  correlation  with  observations  can  be  achieved  with  skewed 
distributions  (either  single-family  or  double-family). 

d.  Two  approaches  are  possible  for  the  derivation  of  the  required  infor¬ 
mation  an  the  strength  of  fabricated  structures.  The  first  involves  the 
separate  assessment  of  the  basic  material  properties  and  of  the  affects 
of  fabrication,  the  hvo  being  subsequently  combined  to  give  the 
resulting  distribution  (the  computer  program  of  Appendix  II  provides 
this  facility).  Tho  second  approach  involves  only  the  statistical 
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analysis  of  largo  numbers  of  identical  components  to  assess  the 
resultant  strength  variation  directly,  without  attempting  to  ascertain 
the  contributions  due  to  the  separate  causes. 

Since  so  much  material  data  exists,  the  first  approach  recommends 
itself,  but  a  deliberate  effort  is  required  to  determine  the  effects  of 
the  various  fabrication  and  assembly  operations  in  statistical  terms. 

12.4  Error  Functions 

a.  The  importance  of  this  function  has  been  illustrated  in  Sections  III 
and  VII,  The  formal  recognition  of  the  probable  discrepancy  between 
the  intended  strength  and  the  actual  achieved  strength  is  perhaps 
more  important  than  the  particular  function  used,  since  the  use  of 
Fayes'  theorem  tends  to  be  self-compensating  once  the  necessary  testing 
is  performed.  The  function  describes  the  discrepancy  (however  caused, 
whether  by  design  errors,  design  tolerances,  deliberate  under-design, 
quality  control  errors  or  fabrication  and  assembly  errors)  in  terms  of 
the  distribution  of  the  probable  mean  strength  of  the  design. 

b.  Section  Vtl  describes  four  types  of  function  suitable  for  the  definition 
of  the  probable  discrepancy.  While  the  use  of  some  standard  function 
is  possible,  it  does  not  permit  the  recognition  of  the  experience  of  a 
particular  constructor  with  his  own  policies  and  practices.  Comments 
on  the  four  types  of  function  ore: 

1)  The  Jablecki  function,  as  used  in  reference  i,  uses  static 
test  data  from  fho  1940  period  (reference  3);  it  is  implicitly 
assumed  that  the  ratio  of  test  strength  to  design  strength  des¬ 
cribes  the  ratio  of  mean  strength  to  intended  mean  strength, 
but  it  is  equally  apparent  that  no  account  is  taken  of  the 
probability  that  the  article  was  weaker  or  st  onger  than 
average. 

2)  Freudenthal,  in  reference  4,  attempted  to  update  the  Jablecki 
data.  The  relevance  of  the  data  used  is  not  altogether  clear; 
the  comment  is  made  that  the  results  ore  representative  of  current 
practice,  yet  doto  are  included  for  aircraft  of  the  pre-1950  period. 
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3)  Both  of  the  above  functions  are  most  easily  used  by 
basing  the  constants  on  a  curve-fit  at  two  selected 
points.  The  same  concept  can  be  employed  using  the 
Gumbel  distribution  of  minima  instead  of  the  Bouton- 
Jabiecki  equation  (linear  log-log  relationship)  or 

the  Freudenthal  exponential  function.  Any  other  suit¬ 
able  distribution  can  also  be  employed. 

4)  The  fourth  type  of  distribution  is  the  double-family  dis¬ 
tribution  described  in  Section  VII.  This  technique  may 
be  the  most  suitable  for  fitting  past  experience,  or  for 
permitting  recognition  of  the  additional  risk  of  design 
error  when  a  radically  new  type  of  construction  is  being 
employed  before  the  required  analytical  tools  have  been 
fully  developed. 

c.  As  will  0  shown  later,  using  the  single-rami ly  Gumbel  distribution 
of  varying  dispersion,  the  degree  of  dispersion  (coefficient  of 
variation)  has  relatively  little  influence  once  the  test  results  have 
been  Incorporated.  A  relatively  low  risk  would  probably  be  intro¬ 
duced  by  the  adoption  of  a  standard  error  function. 

12.5  Presentation  of  Standard  Do  to 
a,  loads  Data: 

1)  Standard  load  spectra  expressed  in  terms  of  some  design 

valoo  (such  as  N-  )  and  of  a  given  shape  con  be  pre- 
Max 

sented  in  tabular  form  as  in  reference  5. 

2)  Mission  profile  and  utilization  data  may  be  standardized 
for  porticulor  aircraft  Or  mission  types,  cuir  will  probably  be 
best  defined  for  each  system  as  port  of  the  specification. 

3)  Data  determining  the  combinations  of  mechanical  and  thermal 
conditions,  the  combinations  of  pilot  and  auto-control  action, 
and  the  combinations  of  external  (gust,  soy)  and  internal  (sub¬ 
system)  effects  cannot  be  standardized  and  must  be  derived  in 
probabilistic  terms  for  each  specific  design.  In  many  coses, 
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this  will  not  be  possible  during  the  design  phase;  some 
standard  arbitrary  distribution  of  effects  may  be  appro¬ 
priate  in  this  phase  for  describing  the  probabilities  of 
engine  failure,  auto-stabilizer  runaway,  cabin  pressure 
malfunction,  etc. 

Strength  Data: 

1)  For  each  type  of  basic  material,  the  present  system  of 
discrete  design  allowables  must  be  retained  for  association 
with  the  design  loads  to  permit  the  physical  sizing  of  the 
structure.  The  values  need  not  be  the  present  "A"  or  "B" 
values  per  se,  but  the  retention  of  these  is  obviously  ad¬ 
visable. 

2)  For  the  reliability  calculations,  the  mean  and  standard 
deviation  (or  coefficient  of  variation)  is  required.  These 
data  are  not  generally  as  readily  available. 

3)  In  addition  to  material  data,  the  statistical  effects  of 
fabrication  and  assembly  are  required.  Typical  values 

for  the  various  processes  (rolling,  stretch- forming,  machin¬ 
ing,  etc.),  for  various  jointing  methods  (riveting,  bolting, 
welding,  bonding)  and  for  the  actual  assembly  process 
(fitting  stresses)  will  be  required,  and  can  be  presented  in 
tabular  form. 

Design  and  Test  Factors: 

1)  For  any  given  set  of  the  other  parameters,  it  is  possible  to 
derive  relationships  between  the  design  factor,  the  test 
factor  and  the  reliability  indicated  by  the  test  result.  Two 
basic  assumptions  will  simplify  the  presentation  in  different 
ways.  The  first  requires  the  adoption  of  a  constant  design 
factor  (soy  1 .5),  but  varies  the  test  factor  to  the  voiue  re¬ 
quired  to  '‘demonstrafeM  the  required  reliability.  Typical 
relationships  are  shown  in  Section  IX. 
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2)  The  second  alternative,  which  may  be  simpler  in  form 
although  less  versatile,  is  to  assume  the  design  and 
test  factors  to  be  equal.  This  retains  the  concept  in 

the  current  system,  but  must  not  be  interpreted  os  having 
the  same  meaning. 

3)  Since  the  "demonstrated”  reliability  is  a  function  of  the 
load  spectrum,  the  error  function,  the  strength  distribution 
and  the  number  of  tests,  it  is  obvious  thot  a  complex  set 
of  charts  must  result.  For  the  particular  choice  of 

o  design  factor  =  test  facto* 

o  one  survival  test 

o  single-family  Gumbel  distribution  of 

maximum  load  per  aircraft  lifetime  (mean 
ot  IDO) 

o  single-family  Gumbel  distribution  of 
minimum  strength  (mean  ot  HX» 

o  single-family  Gumbo!  distribution  nf 
error  {mean  at  100),  where  error  is  the 
ratio  of  achieved  mean  strength  to  in¬ 
tended  mean  strength 

the  curves  shown  in  figures  75,  76,  end  77  the w  the 
manner  in  which  the  factor  can  be  c  Ho  sen  to  realise  o 
defined  reliability. 

4)  Figure  ?S  shows  the  reliabilities  (8)  corresponding  to 
variations  in  ths  load  dispersion  (L  'B  coefficient  of 
variation  of  the  distribution  of  maximum  kW)  and  in 
the  design  and  test  factor.  Separate  carpet  plots  ore 
shown  for  three  levels  of  error  dispersion  (£y),  but  show 
little  variation  with  Ey .  AH  three  corpets  are  tor  strength 
coefficient  of  variation  (Sy)  of  0.04,  and  for  one  survival 
test,  A  series  of  plots  of  this  type  can  be  derived  for  each 
strength  distribution  (for  each  material  type  and  construction 
type). 
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r  -  DESIGN  FACTOR  *  TEST  FACTOR  {ONE  SURVIVAL  TEST) 

Ly  *  COEFFICIENT  OF  VARIATION  OF  MAXIMUM  LOAD 


FIGURE 


DESIGN/TEST  FACTOR  TO  GIVE  CHOSEN  RELIABILITY,  FOR 
STRENGTH  DISPERSION,  VARYING  LOAD  AND  ERROR  OISPI 
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Ly  -  COEFFICIENT  OF  VARiATION  OF  MAXIMUM  LOAD 


Ey  *  ERROR  COEFFICIENT  OF  VARIATION 


F  *' DESIGN  FACTOR  •  TEST  FACTOR  (ONE  SURVIVAL  TEST) 

$v  *  STRENGTH  COEFFICIENT  OF  VARIATION 

FIGURE  76.  OBlGN/TESt  FACTOR  TO  GIVE  CHOSEN  RELIABILITY,  FOR  CONSTANT 
LOAD  DISPERSION,  VARYING  STRENGTH  AND  ERROR  DIVERSIONS 
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F  *  DESIGN  FACTOR  *  TEST  FACTOR  (ONE  SURVIVAL  TEST) 
Ly  *  COEFFICIENT  OF  VARIATION  OF  MAXIMUM  LOAD 


75  OEStGN/TEST  FACTOR  TO  GIVE  CHOSEN  RELIABILITY,  FOR  CONSTANT 
STRENGTH  DISPERSION,  VARYING  LOAD  AND  ERROR  DISPERSIONS 
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P  =  PROBABILITY  OF  SURVIVING  ONE  TEST 
F  -  DESIGN  FACTOR  TEST  FACTOR  (ONE  SURVIVAL  TEST) 


Ey  =  ERROR  COEFFICIENT  OF  VARIATION 

0.04  0.08  0.12 


0.12 


Sv  *  STRENGTH  COEFFICIENT  OF  VARIATION  ALL  LQAD  D1SpERSIONS 


5)  Figure  76  is  o  similar  series  of  plots,  but  the  van  ion 

in  each  carpet  is  with  the  strength  coefficient  of  varia¬ 
tion  (Sy).  The  complete  set  is  for  a  given  coefficient  of 
variation  of  maximum  load.  (Ly  =  0.04).  The  influence 
of  the  error  coefficient  of  variation  is  again  slight.  A  set 
of  this  type  can  be  derived  for  each  of  the  standard  load 
spectra,  and  used  where  necessary  to  guide  the  choice  of 
materia!  or  construction  method  suitable  for  the  attainment 
of  the  required  reliability. 

6)  Figure  77  is  for  a  constant  error  variation  (Ey  -  0.08).  Each 
carpet  shows  the  combinations  of  load  coefficient  of  variation 
(Ly)  and  design/test  factor  and  is  for  a  separate  strength  co¬ 
efficient  of  variation  (Sy).  This  form  of  presentation  is  probably 
the  most  useful  in  the  earlier  design  stages,  when  the  design 
iteration  process  is  being  applied  to  determine  the  layout  and 
member  sizes.  The  example  (figure  77  )  illustrates  the  im¬ 
portance  of  the  strength  variation,  implying  for  example,  that 

a  reliability  of  0.9999  cannot  be  achieved  with  strength  varia¬ 
tions  exceeding  about  0.06,  unless  very  high  design/test  factors 
are  used. 

7)  Figure  78  ‘llustrates  the  associated  probabilities  of  surviving  the 
survival  test.  This  quantity  is  independent  of  the  load  spectrum, 
and  for  the  chosen  equality  of  design  and  test  factor,  shows 
surprisingly  little  variation  with  the  strength  or  error  le  els. 
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SECTION  XIII 

STEPS  TOWARDS  IMPLEMENTATION  OF  THE  PROPOSED  SYSTEM 


13.1  Introduction 

a.  The  proposed  system  of  probabilistic  criteria,  aimed  at  providing 
the  desired  degree  of  static  strength  reliability,  offers  several 
advantages  over  the  present  deterministic  system.  Nevcthele'S, 
it  is  necessary  for  the  essential  features  of  the  system  to  ire  intro¬ 
duced  in  a  m-'-.ner  which  assures  continuity  with  the  p>;  ,ent  system. 

A  two-stage  process  is  suggested  in  this  Section.  Initially,  there 
will  be  insufficient  data  available  to  implement  the  total  aim  of 
establishing  a  single  reliability  figure  covering  the  entire  life  of 
the  fleet;  however,  even  restricting  the  calculations  to  those  flight 
conditions  for  which  data  is  available  will  serve  several  useful 
purposes . 

b.  There  is  an  inherent  resistance  to  new  methods,  especially  when 
the  existing  techniques  appear  to  be  adequate.  Familiarity  (with 
the  old)  breeds  contempt  (for  the  new).  This  is  particularly  true 
in  this  context,  since  the  proposed  method  requires  a  radically 
different  interpretation  of  testing.  Furthermore,  a  number  of 
decisions  will  be  required  which  must  be  based  on  the  correct 
understanding  of  the  probabilistic  processes;  since  this  is  an  un¬ 
familiar  subject  to  many  of  those  who  will  be  responsible  for  the 
decisions,  it  is  vital  that  the  physical,  rather  than  the  mathematical 
interpretation  of  each  step  in  the  chain  should  be  kept  clear. 

c .  The  use  of  the  method  as  a  means  of  comparing  the  relative  risk  rotes 
of  various  designs,  of  various  flight  conditions  and  of  various  struc¬ 
tural  locations  offers  on  opportunity  to  achieve  familiarity  with,  and 
confidence  in  the  method.  It  will  also  encourage  the  acquisition  of 
the  data  required  for  the  eventual  implementation  of  the  complete 
system. 
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13.2  Initial  Implementation 


The  establishment  of  absolute  reliability  values  requires  assurance 
that  every  possible  cause  of  failure  is  considered.  As  this  cannot 
be  guaranteed,  it  is  proposed  that  the  method  be  used  to  establish 
the  separate  probabilities  of  failure  for: 

o  different  structural  designs  under  the  same  loading 
conditions,  in  order  to  indicate  the  optimum  means 
of  securing  the  highest  reliability 
o  the  same  structural  location  for  different  loading 
conditions  (maneuvers,  gusts,  landing,  etc.),  in 
order  to  assess  the  relative  risks  associated  with 
different  flight  cases;  it  is  an  inefficient  design 
which  has  a  high  survival  rate  under  gust  loads, 
but  a  high  risk  of  failure  during  landing 
o  various  structural  locations  under  the  same  leading 
conditions;  this  will  provide  a  means  of  early  assess¬ 
ment  of  areas  of  the  structure  which  will  be  a  poten¬ 
tial  source  of  trouble. 

In  this  context,  it  will  be  possible  to  study  the  influence  of  sub¬ 
system  failures  in  meaningful  terms,  so  that  the  overall  optimum 
can  be  established  for  the  relative  penalties  associated  with  the 
addition  of  redundant  circuits,  or  with  the  addition  of  structural 
weight  to  withstand  the  loads  resulting  from  a  less  reliable  sub¬ 
system.  Information  from  such  studies  will  be  applicable  to  the 
necessary  decisions  which  involve  both  structural  and  non-structuroi 
areas, 

Interfaces  between  structural  design  and  structural  test  decisions 
will  be  studied,  sine*  the  method  provides  information  enabling 
conscious  trade-offs  between  test  load  levels  and  the  probability 
of  destroying  the  test  specimen.  The  necessity  for  testing  to  a 
particular  lood  level  can  be  s'udted  in  terms  of  the  reliability 
level  "demonstrated".  The  necessity  for  redesign  can  be  inter¬ 
preted  realistically  in  the  some  terms. 
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d.  Judies  of  this  type  should  be  performed  on  a  number  of 
existing  operational  aircraft,  us  well  as  on  a  number  of 
new  designs.  This  will  provide  an  insight  into  the  relative 
importance  of  various  parameters,  as  well  as  indicating  the 
implied  reliabilities  of  existing  aircroft  for  the  conditions 
studied. 

e.  It  is  suggested  that  these  initial  studies  should  be  based  on 
the  same  limit  design  factor  as  was  used  in  the  deterministic 
criteria  system.  The  interpretation  of  actual  test  results  will 
be  in  terms  of  the  reliability  indicated  by  the  test  results. 

This  will  provide  the  desirable  continuity  wHh  existing  methods. 

f.  During  these  initial  stages,  it  is  imperative  that  every  induce¬ 
ment  be  given  to  the  collection  and  analysis  of  data  required  by 
the  fuli  method.  This  must  include; 

o  load  spectra  for  different  conditions 

o  ptobabilities  of  different  speed-height-weight 

conditions 

o  strength  distribution  data  for  basic  materials 

o  strength  distribution  data  for  fabricated  components 
using  a  variety  of  fabrication  ond  assembly  methods 

o  achieved  strength  versus  intended  strength  data  to 
verify  the  actual  discrepancy  levels 

13.3  Final  Implementation 

a.  It  will  not  be  passible  to  achieve  a  completely  probabilistic 
svstem  with  any  reel  meaning  until  a  greot  deal  more  statistical 
data  have  been  derived.  However  desirable  o  single  reliability 
value  might  appear,  the  judgment  as  to  what  is  acceptable  will 
remain  arbitrary.  Who  can  decide  logically  whether  0.99996  is 
acceptable  but  0.99994  is  not? 

b.  Because  of  this  dilemma,  it  is  probable  that  the  relative  risk 
assessment  technique  will  prove  to  bo  worth  retaining  even  when 
all  of  the  necessary  data  is  available.  This  provides  not  only  a 
means  of  indicating  potential  sources  of  weakness;  but  also  a  tool 
by  which  the  intended  utilization  con  be  modified  in  such  a  way 
os  t*>  moke  the  best  use  of  o  given  airframe. 
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13.4  Flight  Testing 

One  further  area  which  would  repay  study  during  the  gradual  implementation 
of  the  system  is  the  relatively  high  risk  associated  with  deliberate  flight  test¬ 
ing  to  improbable  corners  of  the  flight  envelope.  The  probabilistic  load 
spectrum  for  such  aircraft  remains  at  a  level  of  1 .0  up  to  the  maximum  in¬ 
tended  load,  which  changes  the  failure  probability  from  that  predicted  for 
the  operational  aircraft.  Studies  of  this  feature  would  probably  enable  a 
more  cost-effective  structural  flight  test  program  to  be  devised  which  is  still 
capable  of  demonstrating  ait  necessary  conditions  at  a  lower  risk  of  loss. 


13.5  Overload  Capacity 

a.  Some  part  of  the  present  factor  of  sofety  has  long  been  recognized 

as  providing  o  margin  of  strength  to  cater  for  occasional  exceedences 
of  the  placarded  limitations.  The  real  overload  capacity  of  an  air¬ 
frame  is,  however,  for  from  consistent,  especially  as  the  structural 
optimization  1$  based  on  the  factored  limit  load  system.  A  frequent 
problem  is  the  solution  of  the  question:  if  a  factor  of  saf'ty  of  1 .5 
exists  at  load  level  P,  ot  what  load  level  does  the  factor  of  safety 
become. 1.0  (or  1,2,  or  1.3)? 


Figure  79  illustrates  the  random  nature  of  the  overload  capacity 
of  the  total  structure.  Suppose  the  external  load,  P,  ot  some  struc¬ 


tural  location  to  vary  linearly  with  the  basic  parameter  (say  N^), 
op.d..*o  po»  through  the  origin  .  os  shown  by  curve  A.  The  un factored 
limit  v dues  ore  N^  and  P^  with  o  factor  of  sofety  of  1  the  design 
load  is  Py  **  !»$  P^,  so  that  the  permissible  with  o  factor  of  sofety 
of  1.0  is  simply  Ny  *  1.5  N^. 


Now  consider  the  existence  of  a  superimposed  loading  independent  of 
N^  (this  might  be  a  locd  system  or  on  internal  pressure,  for 

example).  Curve  8  results  »f  this  loading  odds  to  the  original  loading. 
For  the  some  limit  load  factor,  N. ,  the  unfavored  load  is  now  P. 


and  the  factored  load  is  ?y 
1.0  the  permissible  is® 
the  first  value,  Ny. 


(factor  of  1,5),  For  o  factor  of  sofety  of 

f*ow  Ny  which  is  clearly  greater  then 

8 
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FIGURE  79  DESIGN  LIMIT  AND  FACTORED  LOAD  AND  OVERLOAD  CAPACITY 


2c 


If  the  superimposed  load  system  relieves  the  varying  load,  as 

in  curve  C,  the  loads  P.  and  Py  result,  the  implied  overload 

C  C 

capacity  being  Ny  ,  which  is  less  than  Nn. 

C  U 

Hence,  the  overload  capacity,  Ny,  of  a  given  part  is  dependent 
on  both  the  rate  of  change  of  load  with  loading  parameter  (N^), 
and  on  the  value  at  zero  load  parameter.  If  the  intercept  represents 
a  relief,  the  overload  capacity  will  be  less  than  the  nominal  value, 
but  if  it  adds  to  the  varying  load,  a  greater  overload  capacity  will 
exist . 

Figure  79  represents  the  simplest  of  all  conditions,  a  linear  system. 
The  quantity  which  reflects  the  overload  capacity  will  be  a  local 
internal  load;  in  general,  this  will  not  be  a  linear  function  of  the 
external  load,  end  the  external  load  will  not  be  a  linear  function  of 
any  parameter  which  can  be  used  to  define  the  operational  limitations. 

It  can  be  stated  that  the  actual  overload  capacity  of  a  given  oirframe 
varies  from  one  location  to  another  in  what  is  virtually  a  random  manner. 

The  relationship  between  the  limit  and  omega  (overload)  design  con¬ 
ditions  to  be  used  U  vague.  It  must  depend  on  the  utilization  of  the 
particular  aircraft,  and  on  what  is  regarded  as  a  judicious  risk  of 
failure.  Studies  of  existing  aircraft  should  be  made  to  assess  the 
actual  patterns  of  exceedence  of  limit  condition  and  the  actual  failure 
rotes.  From  such  studies,  it  will  be  possible  to  develop  trends  which 
will  enable  initial  criteria  to  bo  established  which  represent  continuity 
with  present  circumstances. 


27 


SECTION  XIV 

SPECIFICATIONS  AND  HANDBOOKS 


14. 1  General 

o.  The  purpose  of  this  section  is  to  identify  changes  required  to 

MIL-A-8860  through  8871,  MIL-F-8785  and  app'  iate  AFSC- 
DH  series  handbooks  to  implement  the  new  desi^.i  method. 

Implementation  of  the  new  design  method  as  a  replacement  for 
presently  acceptable  procedures  is  not  possible  at  this  time  with 
the  scant  amount  of  appropriate  statistical  information  which 
appears  to  be  available. 

b.  Several  of  the  previous  sections  of  this  report  have  reported  the 
availability  of  statistical  information  and  illustrated  how  it  might 
be  used  to  develop  structural  design  conditions.  In  addition,  it 
is  very  possible  that  much  more  statistical  data  is  available  for 
use  in  the  new  method  than  has  been  uncovered  in  this  brief  study. 
Surely,  many  aircraft  manufacturers  have  in  their  archives  data 
which  is  not  generally  available  concerning  aircraft  they  have 
designed  and  built  and  the  Air  Force  files  undoubtedly  include 
much  data  which  was  not  available  or  not  necessafy  for  use  In 
this  study.  For  example.  Reference  1  implies  that  F-100 
statistics  concerning  vertical  toil  loads  in  operational  usage  ore 
available.  However,  in  this  study  no  vertical  tail  load  statistics 
were  uncovered. 

14.2  Mll-A-8860  Series  Review 

a.  To  start  the  implementation  of  the  new  design  method,  it  is 

proposed  that  appropriate  statements  be  placed  in  the  Mlt-A-8860 
series  (reference  2)  and  in  the  AFSC-DH  series  (reference  6} 
to  allow  the  use  of  statistical  methods  at  an  option.  Then  any 
requirement  for  which  adequate  appropriate  statistical  data  are 
available  can  bo  met  through  the  use  of  those  statistics.  Data 
and  methods  to  be  used  would,  of  course,  be  subject  to  the  ap> 
prove!  of  the  procuring  activity. 
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b.  In  Tables  XX  through  XXIX  the  latest  available  revisions 
to  the  MIL-A-8860  series  ore  reviewed  as  to  applicability  of  the 
new  design  method  at  present  and  in  the  future,  and  data  avail¬ 
ability  to  meet  each  requirement  where  the  new  method  is  applicable. 
Comments  are  included  concerning  changes  required  to  the  subject 
paragraphs  to  implement  the  new  system. 

14.3  Proposed  Changes  to  MIL-A-8860  Series 

a.  In  this  section,  actual  wording  changes  to  the  Mlt-A-8860  Series 
ore  suggested  which  allow  the  use  of  the  new  design  method  as  an 
option .  The  approach  used  results  in  a  near  minimum  number  of 
changes  and  requires  the  use  of  AFFDL-TR-67-1Q7  and  AFFDt-TR-  71 
-178  os  guides  to  incrementing  the  system, 

b.  MH-A-00886QA  (USAF)  31  March  1971 

2.2  Add; 

“ AFFDL-TR-67- 107  Quantitative  Structural  Design  Criteria 
By  Statistical  Methods 

AFFOL-TR-71 -1 78,  Impiementotion  Studies  for  a  Reliability - 
Based  Static  Strength  Criteria  System'4 

3.  Add; 

“Establishment  of  Criteria.  It  is  intended  that  structural 
criteria  be  established  an  a  rational  basis.  Criteria  delin¬ 
eated  in  this  specification  and  the  other  specifications  in 
the  Mlt-A-8860  series  shall  be  used  unless  other  criteria 
are  determined  to  be  more  rational  or  unless  the  criteria  are 
found  to  be  inapplicable  because  of  the  peculiarities  of  the 
aircraft  under  consideration.  New  criteria  or  methods  which 
are  proposed  by  the  Contractor  shall  be  rational  and  shall 
bo  submitted  to  the  U$Af  for  approval  nn’or  to  use  in  struc¬ 
tural  design  computations.  Where  sufficient  statistical 
information  ore  available,  consideration  shod  be  given  to  use 
of  the  methods  of  AFFDL-TR-67-107  ond  AFFDL-18-7’1-1 78 
to  establish  factored  limit  ond  overload  (Omega)  design 
conditions  commensurate  with  prescribed  structural  reliability 
goals. 


TABLE  XX 
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3.4.1  Stotisti col  Methods  Where  approved  statistical  methods 
are  used,  separate  limit  and  overload  (omega)  loading  con- 
dit  ions  and  separate  limit  and  overload  (omega)  factors  of 
sofety  may  be  derived  using  the  methods  of  A.FFDL-TR-67- 
107  and  AFFDUTR-71-1 78. 

3.6  Insert  the  following  after  the  first  sentence- 

“Where  separate  limit  and  overload  (omega)  conditions 
are  derived,  limit  deformation  shall  be  used  with  limit 
conditions  and  overload  (omega)  deformations  shall  be 
used  with  overload  (omega)  conditions. “ 

3-11  Insert  the  following  after  the  first  sentence: 

Mlimit  loads  and  overload  (omega)  loads  shall  include 
applicable  factors  of  safety  where  statistically  aetermined 
limit  and  overload  (omega)  conditions  ore  used.** 

3.12c  Add  the  following: 

“For  statistically  derived  conditions,  allowable  factor  of 
safety  reductions  shot  I  be  negotiated  with  the  procuring 
activity. “ 

6.2.1  Add  the  following; 

“Far  statistically  derived  loading  conditions,  weights 
may  be  established  probabilistically  in  combination  with 
other  design  parameter* .  Weight*  higher  than  the  spec i  - 
fled  maxima  shall  be  considered  in  statistically  establish¬ 
ing  overload  conditions,* 

Add;  “6.2.2.12  Statist  led  Method*  For  statistically  derived 

loading  condition*,  speeds  mny  be  established  prabcbilisticofly 
in  combination  with  other  design  parameters.  Speeds  higher 
than  those  commensurate  with  the  specified  operational  use 
of  the  airplane  shell  be  considered  in  statistically  establishing 
overload  (omego)  conditions.* 

Add:  ”*.5  Structural  Reliability  Cools .  Where  statistical  methods 

of  AFFOi-TR-67-107 or  AFFDl-TR -71  -1  TSare  used,  the  occurrence 
of  limit  on  a  overload  (omega)  lood  levels  ond  minimum  structural 
reliability  goals  shall  be  in  accordance  with  Table  XXX.  ‘ 
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c. 
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3.2  Add; 

“Subject  to  the  approval  of  the  procuring  activity,  the  statistical 
methods  of  AFFDUTR-67-107  and  AFPDl-TR-71-  178  may  be 
used  to  establish  probabilistic  combinations  of  parameters  for  use 
in  the  selection  of  design  conditions.*' 

3.3  Add; 

ac.  Where  sufficient  statistical  information  is  available,  combi¬ 
nations  of  weights  and  load  factors  may  be  established  probabilistically 

3.14  Add; 

"Where  statistical  methods  are  used,  probabilities  of  failure  may  be 
determined  to  establish  levels  of  inoperativeneis  to  be  used  far  design 
conditions," 


3.18  Change  "design  ultimofe"  to  "factored  design  limit  or  facto-ed  design 
overload  (omega)"  in  three  places. 

3.22.2  Add  to  end  of  paragraph; 

**tf  sufficient  statistics  cars  be  established  la  extend  tn*  mission 
analysis  pprooeh  to  omega  toad  extremes,  the  maximum  leads 
derived  from  3,22.2,1.1  alone  mav  be  used  to  gov©™  the  design 
a*  the  airplane,  * 

3.22.2.1.1  Chong*  the  lost  sentence  on  page  l*r  to  read: 

"The  limit  food*  will  be  ndtiphed  by  1.5  to  establish  factored  design 
loads  except  where  siotisticol  methods  ore  used  to  establish  separate 
limit  and  overload  (omega)  loads.  “ 
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d. 

3.1  Add  the  following: 

"Subject  to  the  approvol  of  the  procuring  activity,  the  statistical 
methods  of  AFFDL-TR-67-107  one  AFFDL-TR-71-1 78  may  be  used 
to  establish  probabilistic  combinations  of  the  design  parameters  of 
this  specification." 
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3.2.7  Change  the  last  sentence  to  read: 

"The  analysis  shall  be  performed  in  accordance  with  Mll-T-6053 
except  os  modified  by  approved  statistical  methods.” 


MIL-A-008867A  (USAFJ  31  March  1971 


3.2.2  Change  seventh  sentence  to  read: 

"All  tests  to  design  ultimate  load  (or  to  limit  and  omega  loads 
including  appropriate  test  factors  of  safety  for  statistically  derived 
conditions)  shall  be  completed  prior  to  performing  fail-safe  tests 
and  failing-load  tests  for  any  condition." 


3.4.2  Change  first  sentence  to  read: 

IIT-.i.  i.  -I.  •  l«*  «  I.  I  /I  •  • «  I  l  I*  •  *• 

IQVI9  tv  uii  ut:v4tU  IVVJU  \ftttili  UIlK  VlflC?yu  IUQU5  iHClUUIIig 

appropriate  factors  of  safety  for  statistically  derived  conditions) 
shall  be  .  .  .  . " 


3. 4. 5. 7  Delete  "ultimate-load"  in  the  first  sentence.  Add  the  following 
after  "1 .5  times"  in  the  second  sentence: 

"(or  to  the  appropriate  factor  of  safety  for  statistically  derived 
conditions)" 

3.8  Change  to  read: 

"All  landing  -gear  drop  tests  shall  be  conducted  in  accordance 
with  MIL-T-6053  except  as  altered  by  statistically  derived  landing 
parameters,  where  applicable. " 
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f.  MIL-A-0Q8870A  (USAF)  31  March  197) 

Add  the  following  at  the  end  of  paragraph  3.1: 

"Note:  Subject  to  the  approval  of  the  procuring  activity,  the 
designated  speeds  and  1 . 15  of  this  specification  may  be 
replaced  by  appropriate  statistically  determined  maximum  speeds." 

9  •  MjUA-8871  (USAF)  1  July  1971 

No  specific  changes  to  M1L-A-8871  are  proposed  at  this  time. 

Possible  conflicts  with  implementation  of  statistical  methods  are 
pointed  out  in  tables  XX  through  XXIX. 

14.4  AFSC-DH  Series  Review 

a.  Necessary  changes  to  the  AFSC-DH  series  ir.  order  to 
implement  the  new  procedure  are  quite  minor.  Basically, 
the  changes  involve  redefinition  of  limit-ultimate  load 
concepts  rather  than  use  of  a  1.5  factor  of  safety  in  several 
handbooks  and  the  inclusion  of  definitions  and  reference 
documents  in  DH  1-1 . 

The  proposed  AFSC  DH  1-7,  Aerospace  Materials,  may 
require  some  changes  but  since  it  has  not  been  issued,  it 
was  not  reviewed.  Proposed  changes  for  the  other  documents 
in  the  series  follow. 

b.  AFSC  DH  1-1  (1  December  1970) 

Section  2L,  page  2.  Add  the  following: 

"LOAD,  OMEGA  -  A  low  probability  of  occurrence  over 
load  level  which  replaces  the  ultimate  load  concept  in  the 
application  of  the  statistical  approaches  of  AFFDL-TR-67-107 
and  AFFDL-TR-71- 1 78 

Section  25,  page  1 .  Add  the  following  to  the  definition  of 
SAFETY  FACTOR; 

"In  the  application  of  the  statistical  approaches  of  AFFDL-TR- 
67-117  and  AFFDL-TR-71-178  limit  and  overload  (omega)  con¬ 
ditions  may  have  individual  safety  factors  which  are  s!  Mstically 
determined . " 
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Chapter  4.  Add  the  following  to  the  list  of  references: 

AFFDL-TR-67-107 

AFFDl-TR-71-178 

AFSC  DH  1-.6  (Revised  20  January  1971) 

Design  Note  3BX.  Change  item  3.  to  read  as  follows: 

“3.  Use  an  ultimate  factor  of  safety  of  1.50  except  where 
acceptable  statistical  methods  are  employed  to  develop 
separate  limit  and  overload  (omega)  conditions  and  corre¬ 
sponding  factors  of  safety . " 

AFSC  DH  1-X  (Revised  15  January  1971) 

Design  Note6Al.  Change  item  1.2  to  read  as  follows: 

"1.2  Use  an  ultimate  factor  of  safety  of  1.50  except  where 
acceptable  statistical  methods  are  employed  to  develop 
separate  limit  and  overload  (omega)  condition  and  corre¬ 
sponding  factors  of  safety." 

AFSC  DH  2-1  (Revised  1  October  1970) 

Design  Note  2A1.  Under  paragraph  2 .  BASIC  DESIGN  AND 
TEST  PHILOSOPHY,  replace  the  4th  sentence  with  the  following: 
"Design  the  aircraft  so  that  it  will  not  fail  at  ultimate  loads  (or 
at  limit  or  omega  loads  including  appropriate  factors  of  safety  when 
statistical  methods  of  AFFDL-TR-67-107  and  AFFDL-TR-71-178 
are  used)." 

AFSC  DH  2-X  (15  September  1970) 

Design  Note  1A1.  Change  item  1.2  to  read  as  follows: 

"1.2  Use  an  ultimate  factor  of  safety  of  1 .50  except  where 
acceptable  statistical  methods  are  employed  to  develop 
separate  limit  and  overload  (omega)  conditions  and  corre¬ 
sponding  factors  of  safety." 
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14.5  MIUF-8785B  Review 


a.  This  section  of  the  study  is  concerned  with  establishing  the  need 
and  availability  of  appropriate  data  necessarv  it  the  existing 
design  requirements  of  Ml L-F-8785B  (reference  7)  when  using  the 
new  design  method. 

The  military  specification,  MIL-F-8785B,  contains  the  requirements 
for  flying  qualities  of  United  States  military  piloted  airplanes.  The 
requirements  of  this  specification  should  be  applied  to  assure  that 
no  limitations  on  flight  safety  or  on  the  capability  to  perform  in¬ 
tended  missions  will  result  from  deficiencies  in  flying  qualities. 

The  flying  qualities  of  modern  airplanes  are  the  results  of  in-depth 
design  analyses  using  current  aerodynamic  criteria.  These  flying 
qualities  are  then  evaluated  by  pilots  flying  simulators  or  the 
actual  airplane.  One  of  the  most  acceptable  evaluation  standards 
for  flying  qualities  is  the  Cooper  Rating  System  (reference  8). 

b.  in  MIL-F-8785B  there  exist  three  levels  of  flying  qualities;  i.e., 

Level  1,  Level  2  and  Level  3.  These  levels  are  very  nearly  parallel 
to  the  standards  of  the  Cooper  Rating  System.  The  definition  of  each 
of  the  three  levels  as  specified  in  MIL-F-8785B  is  as  follows: 

Level  1  -  Flying  qualities  clearly  adequate  for  the 
mission  Flight  Phase. 

Level  2  -  Flying  qualities  adequate  to  accomplish  the  mission 
Flight  Phase,  but  some  increase  in  pilot  workload  or 
degradation  in  mission  effectiveness,  or  both,  exists. 

Level  3  -  Flying  qualities  such  that  the  airplane  can  be  con¬ 
trolled  safely,  but  pilot  workload  is  excessive  or 
mission  effectiveness  is  inadequate,  or  both.  Cate¬ 
gory  A  Flight  Phases  can  be  terminated  safely,  and 
Category  B  and  C  Flight  Phases  can  be  completed. 

c.  It  is  not  the  intent  of  this  work  to  regenerate  or  update  the  aerodynamic 
criteria  or  the  flying  qualities  standards.  Rather,  it  is  intended  to 
establish  an  interface  between  these  criteria -standards  and  the  new 
design  method.  This  method  inherently  features  the  statistical  concepts 
of  probability  and  reliability. 
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A  First  round  of  coalescing  these  statistical  concepts  and  the 
flying  qualities  standards  already  exists  in  MIL-F-8785B  and 
in  the  Concorde  flying  qualities  specification  TSS  Standard 
Number  3  (reference  9).  i,;  Hiese  specifications  certain 

degraded  flying  quality  levels  are  linked  with  a  probability 
of  occurrence  of  airplane  failure  states.  No  breakdown  of  the 
failure  states  into  the  various  airplane  components  and  systems 
is  attempted. 

d.  There  are  numerous  aircraft  systems,  such  as  flight  controls, 

powerpiant,  navigation,  landing  gear  and  communication  sys¬ 
tems,  each  of  which  has  different  characteristics  relative  to 
probability  of  failure.  Some  of  these  systems  directly  affect 
the  flying  quality  level  of  the  airplane.  Perhaps  the  most 
directly  related  is  the  flight  control  system. 

In  this  study  the  flight  control  system  of  the  C— 141  MAC  Trans¬ 
port  has  been  chosen  to  illustrate  the  probabilities  of  system  and 
sub-system  failures.  The  C-141  flight  control  system  is  composed 
of  several  subsystems,  the  major  elements  of  which  fall  into  three 
groups;  basic  controls,  trim  controls  and  other  controls.  These 
major  elements  are  further  broken  down  according  to  their  specific 
task  and  they  are  listed  as  follows: 


Basic  Controls: 

1. 

Ai  ieron 

2. 

Rudder 

3. 

Elevator 

Trim  Controls: 

1. 

Roll 

2. 

Yaw 

3. 

Pitch 

Other  Controls: 

1. 

Flap 

2. 

Spoiler 

3. 

Stall  Prevention 

e.  Failure  rate  data  have  been  collected  for  the  C-141  MAC  Trans¬ 
port  fleet  over  a  period  covering  the  entire  flight  life  of  the  airplane, 
which  benan  about  mid  1965.  A  sampling  of  failure  rate  data,  covering 
336,418  flight  hours,  has  been  used  in  this  analysis.  This  data  was 
accumulated  between  September  1968  and  March  1969  and  lists  failures 
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of  each  of  the  above  subsystems.  It  should  be  recognized  that 
this  data  is  but  a  sampling/  that  the  results  represent  trends  and 
are  not  conclusive. 

The  number  of  in-flight  failures  and  in-flight  aborts  due  to  each 
subsystem  of  the  flight  control  system  have  been  extracted  from  a 
voluminous  bank  of  available  data.  The  probability  of  failure  for 
the  various  C-141  sub-systems  is  presented  on  Figure  80.  In¬ 
terestingly  enough  the  trim  control  sub-systems  exhibit  the  lowest 
probability  of  failure  or  the  highest  reliability.  The  boundary 
line  shown  as  Level  2  is  taken  from  section  three  of  MIL-F-8785B. 
The  foctor  used  to  convert  probability  per  flight  to  probability  per 
flight  hour  is  five  (the  nominal  C-141  flight  is  approximately  5  hours 
in  duration). 

Similarly,  the  number  of  flight  aborts  for  each  sub-system  is  shown 
on  Figure  81  .  In  this  sampling  of  probability  data  there  were  no 

in-flight  aborts  attributed  to  the  roll  or  yaw  trim  sub-systems.  It 
should  be  noted  that  the  Level  3  specification  from  MIL-F-8785B  is 
much  more  stringent  than  the  Level  2  standard.  The  scatter  of  the 
data  indicates  that  perhaps  the  specifications  should  be  expanded  to 
cover  separately  each  group  of  sub-systems  such  as  basic  contras, 
trim  controls  and  others.  Up  to  this  point  only  the  flight  con¬ 

trol  system  has  been  discussed.  The  probability  of  in-flight  aborts 
due  to  the  C-141  powerplant  system  for  the  previously  mentioned 
data  sampling  is  72.5  aborts  pet  1000,000  flight  hours. 

f.  Before  on  actual  family  f  specifications  con  be  recommended,  an 
in-depth  study  of  flight  failures  and  in-flight  aborts  is  necessory. 
Typical  classifications  of  airplanes  should  include  transports,  cargo, 
fighters,  tankers,  etc.  Military  transport  and  cargo  fleets  to  be 
analyzed  would  include  the  C-141,  C-5A,  C-130,  KC-135  and 
C-123.  Commercial  fleets  to  be  examined  could  Include  at  least 
the  L- 1 88,  L-10U,  B707,  B727,  B737,  B747,  DC-8,  DC-9, 

DC- 10,  C88Q  and  C99Q. 
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CONTROLS 


FIGURE  81  C- *41  IN-FLIGHT  ABORTS  (336,418  HOURS) 


It  is  possible  that  the  probability  of  failure  analysis  for  each 
classification  of  airplane  could  produce  a  different  set  of 
specifications  for  each.  Even  within  a  classification  the  degree 
of  system  and  subsystem  complexity  can  produce  a  wide  dispersion 
of  failure  probability  data.  In  any  event  it  is  proposed  that  each 
of  these  factors  be  considered  in  the  analysis  to  aid  in  the  develop¬ 
ment  of  a  recommended  set  of  specifications. 


SECTION  XV 

CONCLUSION  AND  RECOMMENDATIONS 


15.1  Conclusions 

a.  The  study  described  in  this  report  has  been  aimed  at  securing  a 
more  complete  understanding  of  the  requirements  for  and 
implications  of  the  static  strength  aspects  ot  the  system  of 
probabilistic  criteria  developed  in  reference  1 .  The  principal 
conclusion  is  that  the  implementation  of  the  complete  system 
would  be  premature,  but  that  a  partia'  application  can  and 
should  be  begun . 

b.  The  concept  of  a  single  numerical  value  for  the  reliability  of 
an  airframe  (or  even  for  one  specific  location  on  that  airframe) 
is  superficially  attractive,  but  any  real  advantage  is  completely 
negated  by  the  problems  associated  with  interpretation  of  the 
number,  Not  only  must  every  porsible  cause  of  loading  be 
established  in  probabilistic  terms,  but  every  factor  affecting  the 
strength  must  also  be  established.  Unless  the  total  picture  is 
assembled  piece  by  piece,  nothing  will  be  known  about  the 
relative  importance  of  the  various  conditions,  and  nothing  will 
be  known  about  ways  of  changing  the  results  by  modifying  the 
operation*  instructions  or  by  redesign  . 

c.  Lock  of  statistical  definitions  of  loading  conditions  is  a  major 
obstacle  to  implementation  of  the  method.  This  is  most  true  of 
asymmetric  flight  coses  and  of  cases  involving  combinations 

of  parameters  (speed,  weight,  load  condition,  load  level,  etc,) 
which  cannot  be  regarded  as  independent . 
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The  reliability  evaluation  depends  on  comparisons  between 
load  and  strength,  both  expressed  by  a  common  parameter.  The 
choice  of  this  parameter  is  complicated  by  interaction  between 
load  systems.  For  example,  if  wing  root  bending  moment  is 
the  measure  of  applied  load,  it  must  also  be  the  measure  of 
strength;  but  the  allowable  bending  moment  may  depend  on 
the  applied  torsion,  shear  and  inlemol  pressure.  Hence  the 
strength  definition  will  generally  be  more  complex  than  implied 
by  reference  1 .  A  normalized  parameter  might  be  used. 

The  need  for  a  single  design  load  remains,  as  does  the  concept  of 
design  allowable  strength.  Without  the  ability  to  match  these 
values,  determination  of  structural  dimensions  is  impossible. 

This  is  recognized  in  reference  I  and  confirmed.  However,  the 
design  factors  to  be  used  will  vary  with  the  statistical  properties 
ioivotved. 

The  strength  distribution  must  recognize  the  variations  due  to 
fabrication  and  assembly  processes,  os  w$|t  as  those  of  the 
basic  material .  Ooto  on  these  effects  is  lacking,  and  b 
urgently  needed. 

The  probability  that  the  ochieved  strength  levels  will  not  be 
those  intended  must  be  recognized  by  the  inclusion  of  a  suitable 
“error"  function  in  the  analysis.  This  toy  be  arbitrary  or  based 
on  appropriate  test  experience;  the  choice  is  relatively  insensitive, 
since  the  incorporation  of  test  results  forms  a  partially  seif- 
compensoting  process , 

Testing  changes  its  roeoning;  it  is  not  a  proof  of  strength,  but 
a  means  of  indicating  probable  error  levels.  The  test  factors 
used  may  vory  according  to  the  reliability  level  to  be  "demonstrated" 
design  to  a  high  foctor  followed  by  testing  too  moderate  factor 
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can  imply  the  same  total  risk  as  design  to  a  moderate  factor  with 
testing  to  a  high  factor.  The  risk  of  destroying  the  specimen 
could  enable  an  optimum  overall  cost-effectiveness  to  be  achieved. 

i.  Repeated  testing  (on  independent  specimens)  will  contribute  to 
the  overall  state  of  knowledge .  Both  laboratory  tests  and  actual 
flight  experiences  hove  the  same  meaning  of  demonstration  of 

a  certain  minimum  strength. 

j.  Test  failures  and  tests  surviving  given  loads  have  different  meanings. 
The  former  are  difficult  to  interpret  consistently,  and  a  test  failure 
should  be  regarded  as  a  test  surviving  a  slightly  lower  load. 

k.  Two  sets  of  design  conditions  require  evaluation.  One  is  aimed 
at  ensuring  negligible  risk  of  a  sample  whore  strength  is  less  than 
the  loads  expected  within  the  placard  limits.  The  second  Is  aimed 
at  providing  a  suitable  margin  of  strength  for  moderate  excesses 
over  placard  limit*,  but  a  lower  reliability  will  be  defined  for 
these  "omega”  conditions.  Different  design  factors  and  different 
tost  factors  may  be  used  for  limit  and  omega  conditions  . 

l.  Management  decisions  will  be  required  which  will  be  based  on 
unfamiliar  information .  It  is  important  that  the  physical 
implications  of  the  various  mathematical  operations  are  maintained 
to  ensure  that  such  decisions  are  correctly  guided.  Compromises 
between  reliability  levels,  design  loads  and  factors,  between 
limit  and  omega  conditton*,  between  weight  end  re  liability,  and 
between  design  and  lest  conditions  w‘-il  be  necessary.  Assessment 
of  the  relative  inpartonce  of  structural  and  non -structural  systems 
will  be  required  in  order  to  achieve  the  requisite  total  reliability 
at  minimum  cost . 


$ e 


GGRESSiVE  VARIATION 


UNDESIRABLE  PRECIPICE 


m.  The  operation  of  the  vehicle  must  be  controllable  in  order  that 
the  intended  reliability  levels  are  achieved.  This  will  require 
selection  of  placard  limits  which  reflect  /he  significant 
parameters,  but  which  remain  practical ;  this  area  requires 
considerable  care .  Subsystem  behavior  may  assume  o  greater 
importance  than  hitherto. 

n.  The  initiol  use  of  the  method  for  relative  reliability  studies  will 
lead  to  the  necessary  acquisition  of  familiarity  with  the  techniques 
and  will  provide  guides  to  the  factors  and  reliability  levels  implied 
by  the  present  criteria,'  gradual  re-evaluation  will  lead  to  more 
efficient  structures  by  permitting  identification  of  structural 
location  and  flight  conditions  which  are  potentially  of  greater  risk. 

o  Continuous  updating  ot  the  evaluations  is  required  to  reflect  she 
increased  knowledge  of  each  stop?  of  the  design  and  operation  of 
a  vehicle , 

p.  The  choke  Of  only  two  levels  {limit  and  omega)  of  which  the 
reliability  gscsls  ore  defined  may  not  always  be  ode  quo  to  tc 
give  the  desirable  progressive  reduction  in  reliability.  For 
example,  increasing  the  temperature  between  limit  and 
<smd§&  conditions  could  result  so  the  onset  of  "thermal  buckling*' 
just  above  limit  condition  and  o  sudden  /eduction  in  allowable  toad. 
Further  increases  might  Have  little  change,  leading  too  "reliability 
praeipke*  of  the  type  shewn  in  Figure  $2  .  It  may  1* 

necessary  to  examine  of  I  cost  one  rendition  between  limit  and 
omega  conditions  to  ensure  the  avoidance  of  such  phenomena  . 

Tht  "high  proof  "  cc.iditicw*  of  reference  I  ?  is  on  example  a* 
this  type  of  precaution  . 


15.2  '-'commendations 


a.  Familiarity  with  the  proposed  system  must  be  gained;  if  is 
recommended  that  a  serie  rudies  be  initiated  which  evaluate 
the  relative  reliability  levels  of  specific  aircraft  and  structural 
locations  for  different  loading  ~ases,  and  of  different  locations 
for  the  same  loading  cases.  The  incompleteness  of  available 
data  is  less  important  in  this  process,  since  a  reasonably  constant 
error  will  have  little  influence  on  the  relative  reliabilities. 

b.  During  this  phase,  attempts  must  be  made  to  collect  and  analyze 
data  which  Is  presently  lacking.  This  includes  statistical 
definitions  of  the  load  sy.  u-r, s  and  of  the  strength  of  fabricated 
structures.  Analysis  of  large  samples  of  existing  (but  relatively 
inaccessible)  test  data  will  permit  selection  of  better  error 
functions  than  those  so  far  proposed. 

c.  The  statistical  equations  used  to  represent  distributions  of  loads 
and  strength  should  be  examined  to  ensure  that  the  important 
"tails"  are  not  required.  Skewed  distributions  and  double¬ 
family  distributions  should  be  investigated. 

d.  The  development  of  the  appropriate  terminology  is  vital  to 
the  understanding  of  the  analysis,  to  the  achievement  of  the 
correct  decisions  for  compromises  and  for  the  selection  of  operational 
guides  which  ensure  that  the  intended  reliability  is  achieved. 

This  new  te.minology  must  recognize  fully  the  changed  meaning 
of  testing;  the  term  "ultimate  loud”  should  be  discontinued  and 
replaced  by  "factored  load";  the  factor  may  be  a  design  factor 
or  c  test  factor  and  the  load  may  be  a  limit  load  or  an  overload 
(omega  load). 
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Initial  application  of  the  proposed  method  to  o  new  design  should 
either 

!)  retain  the  1  .5  design  factor  on  limit  loads  and  vary 
the  test  factor  according  to  the  number  of  tests,  the 
strength  and  load  dispersions,  the  error  function  and 
the  desired  reliability,  or 

2)  use  equal  values  of  design  and  test  factor,  the  value  being 
varied  with  the  same  parameters. 

The  reliability  goals  should  be  based  on  those  implied  by  the  present 
criteria,  to  ensure  no  abrupt  change  in  the  structural  integrity  as  the 
new  method  is  incorporated. 

Specifications  and  handbooks  should  be  modified  to  permit  the  use  of 
probabilistic  methods  as  m  option  to  the  present  methods  where 
sufficient  data  exists. 

The  influence  of  subsystem;  on  the  structural  loads  requires  evaluation 
of  the  rates  of  many  different  lypes  of  failure  .  Acquisition  of 
the  necessary  data  should  be  encouraged. 

Interactions  between  static  strength,  fail-safe  strength  (the 
residual  strength  of  a  damaged  structure)  and  fatigue  "strength” 
require  identification.  Studies  of  the  nature  of  these  interactions 
should  be  pursued  to  permit  the  whole  spectrum  of  structural 
reliability  to  be  expressed  in  a  consistent  manner. 
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APPENDIX  t 

A  NOTE  ON  THE  USE  OF  DOUBLE-FAMILY  DISTRIBUTIONS 


Al.l  It  is  frequently  necessary  to  assume  that  all  of  the  observations  in 
a  sample  are  members  of  a  single  homogenous  population  whose  distri¬ 
bution  follows  one  or  other  of  the  many  standard  forms  (normal,  log- 
formal,  Weibull,  Gumbel,  Poisson,  Pearson,  etc.).  Such  an  assump¬ 
tion  will  often  give  a  good  representation  of  the  observed  probabilities 
of  occurrence,  particularly  in  the  neighborhood  of  the  mode  (the  most 
frequent  values).  For  many  purposes,  a  best  fit  in  this  region  is  desirable, 
but  there  are  other  applications  of  statistical  distributions  where  other 
factors  require  emphasis. 

The  structural  reliability  problem  is  such  a  realm.  The  major  difference 
from  the  more  common  reliability  analyses  is  that  the  "mean  time  to 
failure"  is  not  the  desired  measure  of  structural  reliability.  It  is  the 
risk  of  first  failure  that  is  required,  since  the  ultimate  goal  is  the  prevention 
of  all  failures  (in  effect,  there  is  no  acceptable  failure  rate).  The  implica¬ 
tions  c^e  to  throw  much  more  emphasis  on  the  unusually  high  loads  and  the 
unusually  law  strengths,  which  in  turn  demands  that  the  statistical  repre¬ 
sentations  match  the  appropriate  tails  of  the  distributions  rather  than  the 
regions  near  the  mode. 

A1.2  In  practice,  there  is  no  strict  logic  behind  the  assumption  that  all  members 
of  a  sample  set  of  observations  belong  to  a  single  family,  unless  it  can  be 
verified  that  only  one  independent  parameter  is  involved,  and  this  is  seldom 
if  ever  possible.  Furthermore,  the  information  necessary  to  divide  the  data 
into  its  component  families  will  not  generally  be  available.  Empirical 
methods  provide  a  means  by  which  the  essential  analysis  can  be  performed: 
the  aim  is  simply  to  provide  a  mathematical  model  of  the  population  which 
is  adequate  in  the  region  of  most  importance. 
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A1.3  The  use  of  double-family  distributions  Is  not  r.ev*.;-  power-spectral 

analyses  have  habitually  employed  such  methods,  and  the  representa¬ 
tion  of  loads  and  strength  data  by  two  Gaussian  distributions  is  des¬ 
cribed  in  reference  10  .The  suggestion  that  the  maneuver  loads  spectrum 
may  contain  members  of  two  distributions  is  also  mentioned  in  reference  11. 
This  appendix  expands  the  approach  on  a  more  formal  basis  and  suggests 
methods  by  which  an  acceptable  empirical  distribution  may  be  derived. 

No  attempt  need  be  made  to  ascertain  the  reasons  why  two  families  (or 
more)  are  involved. 

The  examples  are  based  on  the  use  of  the  first  asymptotic  theory  of  extremes 
(Gumbel  distribution,  see  references  12,  13,  and  14)  but  the  principles 
are  applicable  to  any  basic  attribution.  Gumbel's  equations  are  simple 
and  permit  the  easy  formation  of  the  required  quantities  within  a  computer 
program. 

A1.4  Let  the  basic  distribution  be  such  that  the  probability  of  a  value  less  than 
X  is  P,  where  P  is  a  function  of  X,  of  the  mean  (X)  and  the  standard 
deviation  (S)  together  with  appropriate  constants.  In  order  *o  determine 
the  values  of  fhe  constants,  one  viable  technique  is  tc  i-  nsform  the  prob¬ 
abilities  (P)  into  a  new  variable,  Y,  by  means  of  a  transcendental  equation 
which  results  in  a  linear  relationship  between  Y  and  X.  A  least  squares 
best  fit  can  then  be  used  to  match  the  fitted  line  to  the  transformed  observed 
probabilities.  The  pattern  of  the  deviations  is  then  used  as  a  guide  to  the 
choice  of  parameters  for  the  two  families  used  to  achieve  the  desired  repre¬ 
sentation  . 

In  the  case  of  the  Gumbel  distribution,  the  basic  equation  is: 


P  exp  (-exp  (-Y)) 

A 1  —  1 

where:  Y  =  A  +  B 

A)  -2 

A  ■=  ir/d  -  1 . 28255  and  B  -  0.57722 

and  the  transcendental  equation  is 

Y  =  -loge  (-logs  P)) 

Al-3 

A1.5  A  series  of  N  observations  (see  Table  XXXI)  is  arranged  in  ascending  order 
of  X,  each  term  being  allotted  a  rank,  m,  which  ranges  from  1  for  the 
lowest  to  N  for  the  highest.  To  avoid  the  mathematical  dilemma  associated 
with  a  probability  of  one,  the  oc tuct I  observed  probabilities  ( ^jq)  are  replaced 
arbitrarily  by  m/NH  in  the  usua!  manner.  These  values  of  m /N+l  are 

transformed  to  observed  values  of  Y,  using  equation  A 1  —3  and  plotted  against 
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TABLE  XXXI 

DATA  FOR  DOUBLE-FAMILY  EXAMPLE 


X 

m 

m 

N+l 

Y 

240 

i 

.0417 

-1.16 

241 

2 

.0833 

-.91 

243 

3 

.1250 

-.73 

243 

4 

.1667 

-.53 

244 

5 

.2083 

-.45 

245 

6 

.2500 

-.33 

245 

7 

.2917 

-.21 

247 

8 

.3333 

-.09 

248 

9 

.3750 

.02 

248 

10 

.4167 

.13 

252 

11 

.4583 

.25 

252 

12 

.5000 

.37 

253 

13 

.5417 

.49 

256 

14 

.5833 

.62 

256 

15 

.6250 

.76 

258 

17 

.7083 

1.07 

259 

18 

.7500 

1.25 

260 

19 

.7917 

1.45 

261 

20 

.8333 

1.70 

264 

21 

.8750 

2.01 

266 

22 

.9167 

2.44 

275 

23  *  N 

.9583 

3.16  | 

X  =  252.7 
S  -  8.80 

v  =  S/x  =  0.0348 
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as  shown  in  figure  83.  The  best  straight  line  can  then  be  determined 
by  an  appropriate  least  squares  error  method  (reference  32  describes  a  suitable 
technique  which  minimizes  both  the  x-  and  y-  errors).  It  will  be  realized 
that  plotting  the  transformed  probabilities  on  linear  paper  is  simply  equiva¬ 
lent  to  plotting  the  real  probabilities  on  the  appropriate  probability  paper, 
and  for  illustrative  purposes,  figure  84  shows  the  same  data  on  normal 
probability  paper. 

A1.6  It  will  be  noticed  that  the  observations  deviate  from  the  fitted  line  in  an 
ordered,  rather  than  a  tandom  manner,  which  suggests  that  the  assumed 
distribution  is  not  valid.  Now  experience  ..idicates  that  each  of  the  single 
basic  distributions  plots  as  a  line  with  single  curvature  (or  of  course,  as  a 
straight  line);  Vt  is  also  apparent  that  the  data  follow  a  reflex  curve  with 
a  point  of  contraflexure .  The  usual  arguments  as  to  the  importance  of  the 
single  highest  observation  will  apply,  of  course,  and  if  so  desired,  this  point 
may  be  omitted  from  the  best  fit  process.  Even  when  this  is  done,  the 
reflex-curve  pattern  remains. 


A1.7  Now  let  the  assumption  be  made  that  the  data  comprise  representatives  of 

two  families.  Let  these  have  means  and  standard  deviations  XA,  X,,,  S. 

A'  B'  A 

and  Sg  respectively.  Also,  let  Rg  of  the  total  population  be  contained  in 
family  B,  so  that  family  A  contains  (1-Rg)  of  the  total.  The  resultant 
probability  of  a  value  I  ess  than  X  can  now  be  expressed  us 


PT  =  0-Wrb  pb 


AI-4 


where  PA  =  exp  (-exp  (-Y A)) 

Pg  =  e^P  (-exp  (Yg))  Al-5 

represent  the  Independent  probabilities  of  a  value  less  than  X  in  the  separate 
distributions,  where 


VA  '  A 


VA 


x  -  x , 


+  B 

+  B 


Al-6 


The  transcendental  equation  to  derive  the  transformed  probability,  Yj,  is 

the.; 

Yj  -  -loge  (-loge  (Pj)) 
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A 1-7 
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FIGURE  83,.  OBSERVED  AND  FITTED  DISTRIBUTIONS  (GUM8EL  PAPER) 
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Figures  85  through  90  show  the  implied  distributions  in  conventionol 
form. 

It  is  convenient  to  use  the  standard  coefficients  A  and  B  for  both  distri¬ 
butions,  rather  than  to  vary  these;  no  significant  degradation  should  occur 
in  practice,  although  it  would  be  more  correct  to  vary  the  coefficients 
according  to  the  amounts  of  data  allotted  to  the  two  separate  distributions. 

A] .8  The  remaining  problem  is  to  determine  rhe  five  basic  parameters  X^,  S^, 

Rg,  Xg  and  Sg.  Automated  trial  and  error  methods  are  feasible,  but 
simpler  methods  can  be  devised  which  are  generally  satisfactory.  These 
depend  on  appropriate  assumptions  as  to  the  location  of  the  mean  of  the; 
B-family  and  the  nature  of  ihe  overlap.  The  observed  data  are  alloted 
to  suitable  intervals  and  simple  rules  formulated  for  allocating  the  entire 
contents  of  a  band  to  family  B  at  the  upper  end  of  the  range,  allocating 
the  entire  contents  to  family  A  at  the  lower  end  of  the  range,  and  for 
arbiirary  division  between  the  families  for  a  few  bands  close  to  the  assumed 
Xg.  Lockheed-Georgia  Company  has  a  program  of  this  type  which  generally 
provides  good  results,  or  which  serve  as  a  starting  point  for  a  limited  im¬ 
provement  by  trial  and  error.  Once  the  observations  are  allotted  to  the 
two  families,  each  can  be  fitted  by  its  best  straight  line  and  the  compound 
distribution  can  be  generated  from  equations  Al-5  and  Al-4, 

Figures  91  through  93  show  a  worked  example,  using  the  data  of 
Table  XXXI.  The  improved  fit  to  the  observations  will  be  seen. 

A1 .9  The  foregoing  discussion  relates  to  the  case  where  the  distribution  *s  skewed 
to  the  upper  level  of  X.  For  the  opposite  skewness,  the  simplest  way  of 
handling  the  Gumbel  equations  is  to  change  the  sign  of  X  in  the  computations 
(the  derivation  of  a  minimum  value  of  +X  is  equivalent  to  the  derivation  of  a 
maximum  value  of  -X). 

A1.10  It  is  also  interesting  to  note  that  some  observed  distributions  can  be  better 
fitted  by  a  compound  attribution  obtained  by  subtracting  family  B  from 
family  A.  Such  an  approach  may  have  validity  in  strength  estimation,  a 
possiblo  physical  explanation  being  that  the  total  population  consists  of 
several  overlapping  distributions  whose  sum  isc(ose  to  u  single-family  dis¬ 
tribution;  quality  control  processes  then  remove  one  particular  sub-family. 
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TRIBUTIONS  OF  THE  TWO  FAMILIES 


70 


PROBABILITY  Of-  VALUE  < 


P 


FIGURE  92.  CUMULATIVE  PROBABILITY  FOR  D0U8LE-FAM1LY  EXAMPLE 

(G UMBEL  PAPER) 
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Whatever  the  explanation,  the  results.  Inasmuch  as  they  provide  a  good 
empirical  curve-fit  can  be  held  to  be  os  justified  as  the  common  assumption 
that  the  population  is  describable  by  a  single  Gaussian  distribution.  An 
example  of  this  negative  second  family  approach  is  shown  in  Figures  94 
and  95. 

Al-11  For  reference  purposes,  Tables  XXXII  and  XXXlllcontoin  values  of  the 

transformed  variable,  Y,  corresponding  to  various  values  of  the  probability 
of  a  lesser  value  (F)  and  of  o  greater  value  (P). 
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—  DOUBLE  FAMILY 
(FAMILY  B  SUBTRACTED) 


TABU  XXXII 


F 

.99990000 
.99950000 
.99  J00UUO 
. 9950UOUU 
.99000000 
.98000000 
.95Q0OUU0 
.90000000 
.80000000 
.70000000 
.  60000000 
.50(0)0000 
.40000000 
.30000000 
.20000000 
. 15000000 
. 10000000 
.09000000 
.08000000 
. 07UUU000 
• 06U00UUU 
.  ObOoOODO 
. 040 OOOOU 
.1)3000000 
•  02U00U00 
.01500000 
.  01  OOoOOO 
.  0o9U000l) 

.  00800000 
«  Oil 700000 
.  OuoUOOOu 
«  0O50U0UU 
. 0040O000 
.00300000 
*00200000 
.00150000 

•  ooioou  ;o 

.  0009000)) 

« 0iii)80o  )u 
. OU070UOO 
. 0OO60G0O 

.  OUt»uv^UO 

.00040000 


ORDINATES  OF  GUMBEL  EXTREME-VALUE  PAPER 


. OuUlOUUO 
•OU05UOOO 
.OOlOoOOO 
.00500000 
•Oioooooo 

.02QUOOUO 
.05000000 
,  1  ononiioo 
.20000000 
.  30<)OOi):iO 
.40900. n)0 
.  50ot)i)o,)0 
.bOO 00000 

•  70000000 

.  80 1 ;  o :  t <  h  t o 

.  BbOOiiooo 
.90o.I)H)O0 
.91000900 
«92000»o0 
.  930i)i)i /()() 

•  94000001) 

.  98o:ioooo 

•  9bOCiOOUO 
.97000.1  id 
.98000.)  .u 
.  9850Oii.iu 
.  9‘>0  oiiou.i 
. 9910OUUO 
. 99200900 
.9‘>30iJ0U0 
. 994O0OOU 
.  9950o;.. ,0 

•  99<>0OO0u 
.9  o/ouooo 
.998000.))) 
.998500  0 

•  9'  •  ,'0  )!).)() 
.99' '](),•).)  ) 
.99')2U00  ) 
.99930  .9)0 
.99'  ''400  )0 
.°9  'bOli.Hl 
.999600  )0 


-2.220 
-2.028 
-1.933 
-1.6o7 
-1.527 
-1.364 
-1.097 
-.834 
-.476 
-.186 
.087 
.  3b7 
.672 


.0OU3OUOO 

.00020090 

.U00150UU 

.0001001)0 

.00009000 

•Oouoaoou 

.00007000 
.00006000 
.000050')0 
.0001)4000 
.  0OO.l3Jt)0 
.0)101)2000 
•  oi.ooisbu 


.  99970 (JL 
•  99980 UO 
.9998500 
.999900') 
.9999100 
. 9099200 
.9999300 
. 9999400 

.90'. >9500 
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APPENDIX  II 


BASIC  EQUATIONS  OF  MODIFIED  COMPUTER  PROGRAM 

A2-1  Introduction 

A  de$cripf»on  of  the  program  used  for  the  present  study  is  given  in  Appendix  III 
and  examples  are  shown  in  Appendix  IV.  Two  reasons  exist  for  the  use  of  a 
program  different  from  that  in  reference  1 .  The  study  necessitated  gaining  a 
full  understanding  of  the  practical  implications  of  each  step  in  the  procedure, 
and  the  program  of  reference  1  possesses  certain  shortcomings  in  the  extent  to 
which  the  intermediate  results  are  presented.  The  second  reason  was  a  desire  to 
determine  the  degree  to  which  a  given  company  could  utilize  statistical  pro¬ 
grams  already  developed;  the  Lockheed -Georg  io  Company  had  an  operational 
program  for  applying  Gumbel  distributions  in  both  jingle  and  double-family 
form  (references  13  and  1 4),  and  the  incorporation  of  these  was  thought 
desirable.  Combining  these  reasons,  if  was  evidently  easier  to  write  a  now 
program  then  to  odd  fa  the  Original  program  of  reference  1,  although  the  latter 
was  used  as  a  basis . 

A2*2  loads  Spectrum 

a,  Twp  alternative  methods  ore  provided  for  the  definition  of  the  loadprobo- 
bilify;  the  required  form  is  in  terms  of  the  probability  of  o  load  exceeding 
x,»  where  x,  is  a  band-edge. 

b>  Analysts  of  operational  data,  by  the  theory  of  extremes,  provides  o 

suitable  data  base  using  n  minimum  of  information.  If  each  observation  is 
the  maximum  load  in  o  given  rec.r  ding  period  (preferably  a  constant 
period,  such  as  1000  hours),  then  the  distribution  of  such  extremes  is 
&*p*cfed  to  follow  an  exponential  law  of  which  Gumbel 's  equation  is  one 
example.  Reference  12  contains  a  full  description  of  the  theory.  This 
distribution  of  maximum  extremes  is  typically  skewed  with  the  tail  towards 
higher  values,  and  its  use  in  the  present  context  is  suggested  in  reference  1 1 
among  other  sources . 
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Now  the  resulting  distribution  defines  the  probability  (p  )  that  the 

XL. 

i 

maximum  load  exp^  *ed  to  occur  per  1000  hours  (or  whatever  period  is  used) 

is  between  x.  and  x.+dx,  in  other  words,  p  is  the  probability  that  any 
ii  XL. 

]  1 

load  level  up  to  (x.  +  -^dx)  will  occur.  The  resultant  probability  (P  ) 


x. 


of  each  load  (x.)  is  then  given  by: 

i 

P  =  V  p 
x.  Zj  rx. 

j  i=N  i 


A2-1 


as  shown  in  figure  96. 

The  distribution  p^  is  defined  in  terms  of  the  five  basic  parameters  of  the 

double-family  distribution  as  described  in  Appendix  1.  The  cumulative 
probability  of  a  value  less  than  x.  is 


P^  s  exp(-exp(-yA))  •  (1  -  Rq)  +  exp(-exp(-y&)  •  Rg 


A2-2 


where 


^A 


x.  -  x . 

«  1.28255  - - -  +  0.57722 


x,  -  x, 

ond  y.  =  1 .28255  - - -+  0.57722 

8  S8 


A2-3 


when  written  in  terms  of  the  means  ana  standard  deviations  of  the  two 
families,  where,  if  the  coefficients  of  variation  are  defined: 


*.  *  x  •  v. 
AAA 


*8  *KB  *  V8 


A2-4 


If  the  intercepts  and  slopes  of  the  best-fit  straight  lines  on  Gumbo  I  paper  are 
known,  equations  A2-3  con  be  rewritten  os 
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X.  -X 
i  int , 


yA  = 


x.  -  x. 
i  int„ 


yB  = 


A2-5 


where 


and 


x.  „  =x.  -  0.57722  3. 
mt.  A  A 

A 

x.  .  =x  -  0.57722  3r 
mtg  B  B 


3a  =  SA./1 .28255  -  xA  •  vA/l  .28255 
3b  =Sb/1. 28255  =  xg  •  Vg/1 . 28255 


A2-6 


The  cumulative  probability  (P  )  of  a  value  less  than  (x.+dx)  is  similarly 

*i+ 

calculated,  and  the  required  probability  (p  )  of  a  value  in  the  band 

*L. 

i 

x.,  x. +dx  is  found  from  the  difference  of  the  two  cumulative  probabilities. 

i  i 


A2.3  Strength  Distribution 

a.  Material  Strength: 

The  basic  properties  of  the  material  strength  distribution  are  again  input  in 
the  form  described  in  Appendix  I,  but  the  theory  of  minimum  extremes  is 
employed,  which  implies  a  distribution  with  the  tail  towards  lower  strengths; 
this  was  found  to  be  representative  of  actual  data  samples  examined  (see 
Appendix  V)  and  emphasizes  the  importance  of  the  exceptionally  weak 
specimens. 

The  cumulative  probability  of  a  value  greater  than  (x,  +  ~  dx)  is 

Px.  “  exP<'e*P(-yA))  '  0  -*B>  *  =*b(-exP(-yB))  *  Rg  A2.7 

i+ 
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where 


X  —  x  (jx 

yA  =  1 .28255^  s— —  h  0.57722 


*3 


1 .28255 


XB'Xi'2dx 
*B 


+  0.57722 


A2.8 


The  cumulotive  probcbility  (P  )  of  o  value  greater  than  (x.  -  ^  dx)  is 

i  - 

similarly  defined,  and  the  required  probability  (p  )  of  a  value  in  the 

s. 

]  ' 
band  x,  t  jdx  is  calculated  as  the  difference  between  the  two  cumulative 

probabilities.  When  the  whole  distribution  is  defined,  its  overall  mean 

(x  )  and  coefficient  of  variation  (v  -  s,/x  )  can  be  found  by  summation  of 
'v  S  .  s  5 

first  and  second  moments  in  the  usual  way. 


b.  Fabrication  Variation: 

The  available  moteriol  strength  data  above  may  need  modification  to 
recognize  0  secondary  variation  due  to  fabrication  or  assembly  processes. 
This  variotian  is  treated  os  if  is  wore  a  definition  of  the  distribution  of  t;*e 
mean  strength  of  the  materiel,  bet  l  be  the  probability  of  o  men#,  strength 

I  *j 

in  the  interval  x.  tx-dx,  where  the  equations  for  l  ore  pciolle  to  those 

‘  '  /  i 

in  the  previous  paragraph,  The  basic  material  distribution  shape  is  applied 

to  the  fraction  of  the  total  population  which  has  its  mean  c-J  x.  (the  ?  s  dx 

I  2 

range  is  ignored  and  the  sample  assumed  to  occur  at  x,l . 

The  two  families  cro  scaled  so  that  this  fraction  of  the  total  population  is 

!  A2.b 


formed  fr<w  diitributioftt  wi*h  ne ©ft*  o* 


x, 

V  ‘  "a  '  - 

' 


X. 

St  -  v  *  — i 

B.  B  - 

x 

s 


ond  with  the  original  coefficients  of  variation,  ond  v^,  giving  standard 
deviations  of 


*B.  S  *8.  *  VB 
I  i 


A2.10 


The  double  family  distribution  resulting  from  these  values  is  multiplied  by 

the  probability  of  its  occurrence,  namely  L  ,  thus  yielding  the  contribution, 

6p^  ,  to  the  total  probability  of  a  strength  *  x.  $  ^  dx. 
i 


Summation  of  contributions  due  to  all  of  the  mean  strength  voiues  gives  the 
resultant  strength  distribution. 


A2.U 


and  summation  of  the  first  ond  second  moments  enables  the  overall  mean 
strength  and  coefficient  of  variation  to  be  found. 


If  the  secondary  effect  of  the  fabrication  is  not  needed,  this  step  is  simply 
omitted. 


c  .  Where  the  strength  distribution  is  only  retired  os  o  means  of  determining  the 
choree teristic  shape,  and  not  the  absolute  strength  level,  the  units  used 
may  be  chosen  independent);  of  the  units  used  far  defining  the  loads.  The 
necessary  scaling  >»  performed  in  loter  steps. 


A2.4  Intended  Strength 

o.  The  Selected  unfcctored  design  loed ,  which  may  be  either  o  limit  condition 
or  cn  omega  condition,  is  used  os  a  basis  for  determining  the  Intended 
strength  to  result  from  the  structural  sising  procedure.  The  design  factor 
of  safety,  FS,  is  first  opnli  id  to  give  the  factored  design  load: 
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FACLD  =  FS  x  UNFLD 


A2.12 


ond  ony  design  margin  of  safety  then  incorporated  in  the  estimation  of  the 
(factored)  design  load  for  the  present  case 

PDSNLD  =  FACLD  (1  +  MS)  A2.13 

b.  At  this  point,  it  is  necessary  to  consider  the  influence  of  conditions  pre¬ 
viously  exam.ned,  for  the  case  being  currently  analyzed  may  not  be  a 
design  case. 

If  the  previously  critical  design  toad,  DSNLD,  is  greater  thon  the  present 
value,  POSNLD  the  former  is  used  in  all  subsequent  steps,  if  the  new 
value,  PDSNLD,  exceeds  the  previous  value,  it  replaces  DSNLD. 

lh« strength  levels  implied  by  other  structural  constraints  such  as  stiffness  or 
fatigue,  ore  incorporated  in  the  some  way,  on  appropriate  vgluO  of  DSNLD 
being  input, 

c.  Once  the  critical  design  bod,  DSNLO,  is  established ,  it  is  related  to  the 

design  olbwoble  strength,  expressed  os  o  number  of  standard  deviations 

below  the  inferwed  mean  strength  If  the  conventional  'A*  value  is  being 
used,  this  will  be  roughly  2.33  (assuming  »ha  strength  distribution  is  norma!}. 
It  therefore  follows  that 

OSNU,  -  AMSTR  (1  -  5  .  •  VJ  A2.U 

alt  S 

where  AMSTR  is  the  intended  mason  strength 

is  the  coefficient  of  variation  of  the  strength  distribution. 

Hence 

AMSTR  DSNLD /(l  -  S  n  -  VJ  A2JS 

oH  S 

which  defines  the  intended  meon  strength  of  the  total  produc  tion  run  of  the 

part  being  analyzed. 
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d.  No™  the  equations  of  section  A2.3  define  the  distribution  of  a  population 

whose  mean  is  at  x  .  The  actual  mean  is  intended  to  be  at  AMSTR,  and 
s 

the  intended  distribution  of  strength  of  the  individuals  in  the  production 
run  is  obtained  by  repeat snr  these  steps  with  and  x^  values  replaced  by 

XA  *  AMSTR/ xs 

0n<i  x„  '  AMSTR/x 

b  5 

respec  lively . 


A2.5  Intended  Reliability 

a.  With  the  strength  distribution  resulting  from  the  previous  paragraph,  the 
probability  that  the  strength  lies  in  the  interval  x  *  i  dx  is  known  os  p 
A  structure  Having  this  strength  will  fail  if  the  load  exceeds  x.  Si 

(this  Is  not  strictly  true,  04  it  will  not  foil  at  load  x.  if  th*  strength  is  in 
the  upper  half  of  the  interval;  however,  this  necessary  approximation  intro¬ 
duces  negligible  errors  if  the  intervol  sire,  dx,  is  not  too  large).  The 

probability  H*ot  the  load  exceeds  x.  is  already  known  to  bo  P  ,  hence  the 

'  *1. 

I 

probability  of  failure,  which  represents  the  simultaneous  occur -gnee  of  these 
two  events,  is 


M>r  =  P  *P 
f.  x  x 
«  ».  I, 

I  t 


A2.16 


b.  Integration  over  the  whole  range  of  strength  yields  the  total  risk  of  failure 

N 

A2.17 


p^  X  *p 

f  ~  f 


ii 

and  the  reliability  is  the  complement  of  this,  namely 

R  -  I  .  Pp 


A2.I8 


c.  One  simplification  con  be  mode  for  computation:  If  the  strength  distribution 
is  summed  to  give  the  probability,  P  ,  that  the  strength  is  less  thon  x,: 
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A2.I9 


W-  F 


=  2 


'R  Z,  px 

sj  i=l  si 

then  up  to  the  highest  load  level  for  which  P  =1.0,  the  failure  proba- 

XL. 

I 

bility  can  be  expressed  in  one  step  as 

P_  =  PD  -1.0  A2.20 

r .  K 

1  si 

and  integration  can  be  started  at  this  l^vel. 

d.  The  significance  of  the  &Pp  values  and  of  the  cumulative  integration  of  Pp 
ore  of  some  interest.  The  distribution  of  6Pp  indicates  the  density  distri¬ 
bution  of  the  risk  of  failure  and  can  show  whether  a  greater  gain  in  reliability 
could  be  achieved  by  operational  restrictions  or  by  deliberately  modifying 
the  strength  distribution.  A  peak  at  low  x-values  indicates  that  the  very 
weak  specimens  ("certain11  to  fail  because  of  the  high  probability  of  the 
load)  contribute  most  of  the  total  risk.  A  peak  at  high  x-values  indicates 
that  the  rare  high  loads  are  the  major  cause  of  the  total  risk,  in  the 
former  case,  little  gain  would  lesult  from  elimination  of  high  load  levels, 
but  in  the  latter  case  the  benefits  would  be  greater. 

A2.6  Probable  Discrepancy 

a.  The  probability  that  a  discrepancy  may  exist  between  the  intended  strength 

and  the  actual  strengt1  of  the  design  is  next  incorporated.  Algebraically, 

this  is  performed  by  means  of  an  assumed  distribution  of  achieved  mean 

strength,  p  .  The  program  of  Appendix  III  contains  four  alternative 

SM. 

i 

functions  suitable  for  this  purpose. 

b.  Bcuton/Jablecki  Function; 

Reference  1  describes  the  equation  used  to  represent  the  test  data  accumu¬ 
lated  by  Jablecki  from  tests  performed  during  the  1940  decade  (reference 


&  I 
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3) .  The  equation  represents  a  linear  variation  of  the  cumulative  proba¬ 
bility  of  failure  with  the  ratio  of  achieved  load  to  intended  ultimate 
strength,  both  being  plotted  on  logarithmic  paper.  The  program  of  reference 
1  locates  the  upper  end  of  the  line  at  a  point  with  a  cumulative  probability 
of  1 .0  when  the  load  ratio  is  1 .185,  and  in  the  "standard"  case,  locates 
the  lower  end  at  0.01  probability  when  the  load  ratio  is  0.333.  Other 
levels  are  defined  by  varying  the  probability  of  failure  at  this  same  load 
level  of  0.333  (figure  97). 


The  present  program  retains  the  same  general  function,  but  permits  the 
input  of  any  two  points  on  the  straight  line.  The  maximum  cumulative 
probability  is  truncated  at  1 .0.  The  resulting  equation  is  then  used  to 
generate  the  probability  than  the  mean  strength  lies  in  each  of  the 
intervals  x.  ±  ^  dx.  The  equations  used  are  as  follows: 

Let  PFj  be  the  given  probability  of  failure  below  load  PPUj 

andPF^  be  the  given  probability  of  failure  below  load  PPU^ 

then  the  general  probability  of  failure,  PF,  below  load  PPU  is 
given  by 


log10PPU-log]0PPU1  log10PF-log10PF1 

T°9l0^U2TTog10PPUl 

whence  log10  PF  =  RI  jlog^ Q  PPU  -  log^  A 


where 


and  log^Q  A  = 


Io9I0PF2“1o910PF1 

KI-|og-0PPU2-!o9l0PPU, 

_  lo9l0  PF2  1o910  PPU1  ‘  '°910  PF1  lo910  PPU2 


log10  PF2 


,ogioPhi 


A  =  10 


log10  PF2  l°8lOPPUl  '  Io9|Q^1  ‘°91QPPU2 

l°S)otV',°3:o  PFi 


A2.2I 


A2.22 


A2.23 


A2.24 
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hence 


So  PF  “  loEI0  '  A  / 


or 


PF 


(ppuy 

Rl 


Ri 


A2.25 


The  fraction  of  the  total  mean  strength  distribution  lying  in  the  interval 

x.  ±  i  dx  is  then  found  from  the  difference  of  the  cumulative  probabilities 

1  2  11 

x  -  ^  dx  x  +  ^  dx 

at  the  two  band  edges,  substituting  ar>d  ~pSN~LD~  'n  *urn  ^°r 

in  equation  A2.25  with  A  and  Rl  determined  from  equations  A2.24  and 

A2.23,  respectively.  The  resulting  differences  give  the  required  values  of 

p<-  for  each  interval . 

M. 

i 

c.  Freudenthal  Function: 

The  expression  used  is  a  general  form  of  the  equation  on  figure  3  of 
reference  4: 


p  l  PPU( 

ps "  ®*p  r 


Rl 


A2.26 


where  P<.  is  the  probability  of  exceeding  a  mean  strength/design  strength 
ratio,  PPU.  The  corresponding  probability  of  drawer  value  is 


pc  -  i  /  PPU  \ 

PP  ]  'exP  (“~a 


,RI 


A2.27 


As  with  the  Jablecki  function,  the  present  program  enables  the  constants  to 
be  derived  from  two  known  values  of  PF  and  PPU.  At  the  two  given  points, 


PF .  =  1  -  exp 


PPU 


,RI 


PFj  -  1  -  exp  (* 


PPU, 


xRI 


A2.28 


whence 


90 


'PPU, 


vRI 


—  =-IogeO  -PF,) 


'PPU, 


vRl 


-loge  (1  -PF^ 


or 


Rl(log10  PPU]  -  log,Q  A)  =  log,0(-loge(l  -  PF,)) 
RWo9,0  -  log10  A)  -  Iog10(-lo9e(l  -  PF^) 

Solving  for  Rl  yields 

iog]0Hoge(1  -PF^)  -log10(-loge(l  -PF,)) 

Rl  ’  B?t0wVKom,i  - 

and  sobsi  uting  for  A  gives 


A2.29 


A2.30 


A  “  PPU,/(-loge(l  -  PF,)) 


l /HI 


A2.3I 


The  particular  system  routine  for  calculating  logarithms  prevents  zero  or 
unity  being  chosen  as  input  values  of  PF  or  PPU. 

The  procedure  for  evaluating  the  values  oF  is  identical  to  thot 
previously  described.  i 

d.  Gumbel  Function: 

The  process  is  essentially  the  some  as  the  above,  but  the  Gu-nbel  distribu¬ 
tion  function  of  minimum  extremes  is  used  to  fit  a  tine  through  the  two 
input  points  which  ore  defir^i  as  before.  At  tho  two  given  points,  the 
probabilities  of  lower  voiues  ore 


PF,  4  1  -  exp(-exp( -y,)) 
P?2  “  *  “  «xp(-exp(-y2)) 


A2.32 


where 
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A2-33 


PPU  - PPU 


yjaA - S - 

PPU  -  PPU, 

y2  =  A  S  * 


+  B 


where  A  =  3 .28255 
B  =0.57/22 

PPU  is  the  mean  of  the  implied  distribution  and  $  is  the  standard 
deviation . 


ppu.  - m 

Now  A  =  A - =- - B  =  log  {-log  (1  -  PF,)) 

l  5  eel 

PPU  -PPU 

and  A  =  A - = - -  B  =  log  (-tog  (1  -  PfJ) 

l  >  e  e  & 


A2.34 


whence 


PPU 


{A]  *  B)PPU2  -(A^  B)PPUt 

A1  "^2 


A2.35 


and 


$  - 


Aj  ♦  B 

A(PPU  j  -#U) 


A  2. 36 


The  calculation  of  the  distribution  of  mean  strengths,,  p<.  ,  is  then  os 

follows:  i 


5 


The  probability  of  a  value  greater  than  x.  -  ^  dx  will  be 
Pfcj  =  expt-ex^-y^)) 


•  , 

,  ..  » .282$S  itasn  *1*1  °*  l . 

where  y^  4  — -g - IrPU  •  j bfc"  '  0.57222 


A2.37 


»  «. 


j 

and  the  probability  of  o  value  greater  than  x.  ■»  ?$■  dx  will  bo 

PR  2  “  exp(-exp(-y2)) 


A2.38 


where 


-wirh0'57722 


so  thqt  the  required  population  in  the  interval  x.  ±  -=  dx  is 

i  2 


PS  =  PR,  -  PR2 


A2.39 


e.  Double-Family  Gumbel  Distribution: 

The  fourth  error  function  permits  the  use  of  the  more  general  double- 
family  distribution  defined  by  the  means  of  the  two  families,  the  standard 
deviations  {these  ore  assumed  equal  in  this  application,  so  that  the  number 
of  input  parameters -remains  at  four,  os  with  the  ■sther  error  functions),  and 
the  fractions  e#  the  total  distribution  allotted  to  each  family  (the  fraction 
belonging  to  5>he  tower  strength  family  is  actually  input). 

The  equations  for  the  distribution  ore  similar  to  those  described  in  Section 

A2.3,  but  with  tho  probability  (p  )  replaced  by  die  probability  (p,  )  of 

*s.  V 

»  i 

each  mean  strength. 


A  2J 


th,  wit 


a .  The  district. -an  of  mean  strengths,  p<.  »  which  defines  the  assumed  error 

M, 

i 

function,  ii  combined  with  the  basic  indb'idual  strength  distribution,  p  , 

s, 

i 

to  give  the  probable  individual  strength  distribution  in  the  presence  of  the 
wrar.  The  synthesis  process  is  idenficoi  to  that  described  in  section  A2.3.b 
for  the  inccrpccotion  of  the  fabrication  variation.  The  new  distribution, 
p  ,  replaces  the  previous  distribution,  end  its  cumulative  version,  P  , 


is  formed  to  replace  the  previous  vofues  of  P,t 

S. 


b.  The  probability  distribution  of  failure,  6Pp  ,  and  the  cumulative  probability 

i 

of  failure,  Pp  ,  are  re-estimated,  following  the  procedure  of  Section  A2.5, 

i 

enabling  the  reliability  to  be  estimated  for  the  revised  state  of  knowledge 
(with  probable  discrepancy,  but  before  testing). 


A2.8  Incorporation  of  Results  of  First  Test 


a .  The  next  set  of  updates  can  either  predict  the  effects  if  certain  test  results 

are  assumed  to  occur,  or  can  be  used  to  revise  the  estimates  after  actual  test 

results  have  been  obtained.  Before  this  step  is  performed,  it  is  useful  to 

predict  the  chance  that  the  first  test  load  will  be  survived.  The  first  test 

load,  x_  ,  is  defined  as  the  unfactored  load,  UNFLD,  multiplied  by  the 

1 

desired  test  factor,  TF, .  Now,  the  cumulative  probability  that  the  strength 
is  less  than  x.  14  refined  as  P„  . 

1  Rs. 

I 

Hence,  the  probability  thot  the  strength  of  the  first  test  specimen  will 
exceed  the  test  load  is 


I 


where  j  satisfies  tl  <  condition 


X ,  -  r  dx  <  X»  ■i  X .  +  ~  dx 

I,  2  Tl  h  2 


where 

x^  *  UN  FLO  •  TP, 


A  2 -40 


A2-41 


A2-42 


and  Pj.  is  then  the  required  probability  of  surviving  the  first  test. 


The  probability  of  o  second  specimen  surviving  o  tost  to  o  load  given  by 


xJ  *  UNFLD  •  TF2  A2-43 
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is  similarly  calculated  to  be 


where 


x.  -  x  dx  <  x_  <  x.  +-?rdx 

l2  2  T2  l2  2 


so  that  the  probability  of  surviving  both  tests  is 


A  2-44 


A  2 -45 


F 


T 


2 


A -46 


and  the  same  process  is  repeated  for  the  required  number  of  independent 
tests  (of  different  specimens). 


A2-47 


b.  Three  different  types  of  testing  can  be  selected.  The  first  procedure  consists 
of  performing  tests,  each  surviving  the  same  load  level.  The  second 
consists  of  a  series  of  tests  to  failure,  each  of  the  tests  being  to  a  dif¬ 
ferent  load,  XT  .  The  third  possibility  is  a  series  of  N_  survival  tests, 

N  1 

each  test  surviving  a  different  load  level,  X_  ,  The  implicc'  o  have 

N 

been  discussed  in  Section  IX,  and  this  Appendix  will  simply  give  the  equa¬ 
tions.  Bayes'  theorem  is  used  to  modify  the  distribution  of  mean  strengths 
in  a  manner  which  reflects  the  test  results  (assumed  or  actual) .  Reference 
15contoins a  useful  example  of  this  particular  application. 

c.  Nj  tests  surviving  the  same  load,  X^.: 

The  prior  (before  test)  distribution  of  mean  strength  is  already  known  to  be 

p,  .  The  posterior  distribution  Is 
M. 
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A2.48 


'  Pc 


1P$ 


*T/x=x; 


M. 


M. 


1  Pc 


T  /- 
'  x=x. 


M. 


where  (P<.  )~=x  is  the  probability  of  surviving  the  test  load,  X^,  when  the 


XT  ' 


strength  distribution  has  its  mean  at  x.  ±  ^  dx,  p<.  is  the  probability  that 

'  1  v 

1  i 

the  mean  is  at  x.  ±  dx,  and  the  denominator  is  a  normalizing  factor  to 
ensure  that  the  total  posterior  probability  of  all  mean  strengths  remains  at 
unity.  The  values  of  are  calculated  from  the  dispersion  properties  of 


the  basic  strength  distribution  with  the  actual  values  scaled  to  give  a  mean 


at  x. . 


If  N-j.  is  greater  than  one,  the  process  is  repeated,  later,  but  with  the  pos¬ 
terior  distribution  of  the  previous  iteration  used  as  the  prior  distribution  for 
the  subsequent  iteration. 


It  should  be  noted  that  the  resultant  effects  of  this  test  procedure  are  identi¬ 
cal  to  those  used  in  reference  1 .  Equation  A2-48  reduces  to  the  following 
form  if  the  denominator  is  assumed  to  be  unity; 

T 

•  p,  A  2 -49 

M. 

I 

x=x. 


d.  Nj  Tests,  Each  Failing  at  Load  : 

i 

A  similar  process  is  employed,  but  the  probability  of  surviving  the  tost  load 
is  replaced  by  the  probability  that  the  strength  lies  in  the  interval  con¬ 
taining  the  test  load.  Hence,  tho  posterior  distribution  of  nuxm  strung!  vs, 
after  the  first  test,  is 
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id  so  on .  The  dependence  on  the  interval  width  is  evident. 

e.  N.j.  Tests,  Each  Surviving  Load  x^ 

i 

The  procedure  is  similar  to  the  first  procedure  (paragraph  c)  but  o  different 
x-j.  is  used  for  each  posterior  condition . 

f.  Whichever  test  process  is  employed,  the  revised  estimate  of  the  distribution 

of  mean  strengths  is  again  used  to  revise  the  distribution  of  individual 

strengths,  enabling  a  new  definition  of  p^  .  This  is  used  in  turn  to  re- 

x. 

evaluate  the  failure  probabilities  and  the  reliability,  following  the  steps 
described  in  A2.5. 

A2.9  Incorporation  of  Subsequent  Tests 

a.  After  inclusion  of  the  first  testresu.i,  a  revised  estimate  con  be  mode  of 
the  chance  of  surviving  the  second  and  further  hs»H.  The  equations  are 
equivalent  to  those  in  Section  A2.8.o. 

b.  The  means  of  revising  the  probable  distribution  (p^  )  of  mean  strengths 

M, 

l 

Isos  been  described  in  the  previa- •»  section  .  The  appropriate  posterior 

distribution  leads  to  on  updated  distribution  (p^  )  of  individual  s?r«ir.gfhsf 

x. 

i 

end  so  to  a  revised  estimate  of  the  failure  risk  and  of  the  reliability. 
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APPENDIX  III 

COMPUTER  PROGRAM  USED  IN  STUDY 


A3. 1  Introduction 

The  program  used  in  the  study  was  based  on  the  STTREL  program  of 

reference  1.  The  modifications  desired  made  a  new  program  easier  to 

write  than  their  incorporation  into  the  existing  program.  These 

modifications  comprised: 

o)  step-by-step  computation  and  print-out  of  the  vorioos  stages 
of  the  total  procedure 

b)  a  constant  calculation  interval  to  clarify  interpretation 

c)  the  facility  for  superimposing  o  fabrication  variation  on  to 
the  basic  material  strength  distribution 

d)  o  wider  voriety  of  error  functions 

e)  the  facility  for  assessing  foi lure  tests  and  survival  tests  to 
different  test  levels 

f)  the  use  of  double -family  Gumbei  distributions  throughout,  except 
that  load  spec*rym  ordinotes  eon  be  input  in  place  of  this  dis¬ 
tribution. 

A3. 2  Summery  of  Fre^nom 

o)  The  system  comprises  o  moin  program  (STPRj,  eight  subroutines, 
on*  function  end  q  dato  block  sub-program.  A  flow  chart  of 
the  main  program  appears  in  figure  9S.  A  brief  description 
of  the  mo>n  program  and  of  each  of  the  subroutine*  follows  in 
conjunction  with  listings  of  the  source  decks.  The  logic  employed 
was  at  simple  os  possible,  in  the  interests  of  clarity  end  no  attempt 
was  mode  to  minimize  run  limes. 

The  program  is  written  using  FORTRAN  V  for  the  UNIVAC  1 106 
Computer  with  the  EXEC -8  operating  system.  A  CAICOMR  plot 
option  is  available  end  requires  one  magnetic  tope  when  used. 

VVhea  the  plot  is  not  usu-,  the  only  peripherals  required 

ere  the  card  reader  end  printer. 
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SET  STANDARD  DATA 
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(b)  INTENDED  STRENGTH  STAGE 

FIGURE  93,  FLOW  CHART  OP  MAIN  PROGRAM  (CONTINUED! 
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(c)  PROBABLE  DISCREPANCY  STAGE 

FIGURE  98.  FLOW  CHART  OF  MAIN  PROGRAM  (CONTINUED) 
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(d)  TEST  RESULT  INCORPORATION 

FIGURE  98.  FLOW  CHART  OF  MAIN  PROGRAM  (CONCLUDED) 
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TABLE  XXXIV 
STPR  INPUT  LIST 


TABLE  XXXIV  (Continued) 
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JOS 


140  AJ'LCT  0.0  Plot  Option:  0.0  a  No  plots 

1.0  **  Plots  required 


One  note  regarding  output  must  be  made;  the  program  was  also 
operated  on  the  multiple  terminal  remote-access  (DEMAND)  system 
in  use  at  Lockheed-Georgia  Company.  This  system  possesses  two 
output  modes;  WRITE  (6,  XXX)  causes  output  to  be  printed  on-line; 
WRITE  (10,  XXX)  enables  the  output  to  be  internally  stored  for  later 
offline  display.  The  code  NP  is  set  to  6  or  10  according  to  the 
sign  allotted  to  the  cose  number.  The  option  may  be  easily  changed 
to  suit  the  available  output  deviee  codes. 

A3. 3  Input  Data 

This  is  defined  at  this  point  since  use  nf  the  defined  operating  controls, 
etc.,  simplifies  the  program  descriptions  which  follow.  Table  XXXIV 
defines  the  various  items  with  their  locations  in  th_  data  block  and  the 
standard  values  built-in. 


A3. 4  Description  of  Program 

A  listing  of  each  of  the  routines  is  given,  with  notes  describing  the  purpose 
of  the  appropriate  section.  The  basic  equations  are  given  in  Append!*  II, 


If  the  plot  routines  ore  unsuited  to  the  user's  computer  system,  input  item 
140  shr  id  be  ignored  and  the  following  cards  removed: 

65  through  6?,  122,  134  through  138. 

317  through  321,  335  through  338,  350  through  354, 

403  through  410,  419  through  422,  42?  though  43?, 

471  through  477, 

731  through  738, 

802  through  809, 

994  through- 1000,  1017  tKough  3030,  1045  .hrough  1058, 

a)  MAIN  PROGRAM 

The  first  69  lines  of  the  program  control  the  allocation  of  storage, 
the  definitions  of  common  and  data  blocks,  equivalence  statements 
and  other  system  controls. 


5  through  35, 
299  through  309, 
373  through  386, 
446  through  453, 
560, 

777  through  787, 


482  through  493  ,  503  through  514, 
?40  through  743,  75 8  through  766, 
812  through  90i,  905  through  992, 


106 


_  1‘ 

2* 

J* 
n  " 

S  * 


_ Ml*.  _ _  _  _ _ __ 

OloEhSlOfe  NAHM  I?  I  ,T*tai  ,P«S21200) 


_$tPKQ0CJ_ 
STf-x  J002 


COhHQW  />/  PAL  4  3001  , « 1200  )  .  SX  ,  P«S  I  200  )  .  PS  5  I  200 )  .PSH  1 300  I ,  PSh2  1  200  )  STpXOOu  J 
C  *  •  •  «  S  T  P  POLL,  ¥ 

C _ COPHOH  rox  CALCOHP  PL 0  Ts _  STPkOOJS 


«• 

_ It 

4< 

«• 

To. 
1 1  ■ 


C01IH0*.  /COX^ir/  U'LTT  1*02  >  .  tPLTT  l*H)2>  ■  VPLT842Q2) ,  *PcTB<202) 
J_  »HUF*,rt05C  I_26  ,f  1  02.F22  .M.  I  N  I  P .  HL  1  N2P  , _ 

J~  ><Cin>t.hi.i  s  ,  e  b'J>  tioY 

__p  1  Nl'nSt Ob  1 0  U  P  4  S  0 1  _  _ _ _ _ _ _ _ _ 

E  -Vu  l  V  A  L  £  b  C  C  It  64)  F  (i  )  » L  8  41  F  4  l  l  ')”" 

QIH^nS  10s  IF  I2AUI.  1KU2I  10) _ 


$T  P, 4juu* 
SIPKOLij? 
"STPhlJOuU' 
S'PPjOy 
■'SIPHOufu 
$T?„0l)| 1 
TTphouTT 
STPVOO 1 0 
'SlPAUUi  s' 
STpwOUlb 
STphUs  1  »' 
BASIC  STpWoOi; 


12* 
J  3*_ 

l*t» 

I  £>• 

"l** 

1  ?• 


c* « • « 
JL 


aufrsjs  row  PLOT  TITLES 


OlMEwVtOti  64)>*  *4  II  21  ,h00C1  IV)  ,ri2»  ti)  ,F|12||  it  I 
1  12216AI  ^MLIlilP(A)  iMLibiPtSI  iHLl'i.iPiai  iHLIittHH 
o  1  oEns  io»J  r  It  i  ti  2i  .r  1 2 1 1 1 2 1 

DATA  4  F  i  I  I  1  I  1  .  I  «  I  ,  1 1 1  /*  il.Tl.bOED  ST»CB(.Th  - 


4«* 

ay 

21* 

'i2«' 

20* 


i*CAN  STKEnLIm  «  •/ 

__OJ-«5>>tO«  ....  F  »  ist  I  I  )  .F  1201  I  I  1  ,E  I  31  I  I  n  ,m».4  1 1 ) 

C-tlTi  VAlEJiCC  TnOCC  (Vl  .NAME  tm' 


•wnurnr 

SlPMQOlV 

'STpHOdio 


DATA  IFIlSllm»l  .  II  >/«  UtCNPCO  FAlUlKE  »WU8..  BQ  0  |  SC»CP  A  NC  T  .  _KO  STPP002I 

...  ..  ■  '  SV*»M30>>‘ 


ircsr 


1**LT  i  At>  TEST 

*  / 

STpuullJb 

25* 

OAT*  IF121  11  1  ,  lot  ,  12)  /•SLVlSlO  HEAtt 

ST»(>gTm 

■ 

STpBQf.FS 

l  VAk  * 

•  / 

Slp,OU2» 

TU.LA _ ‘iLAlXl /»liPt)Attf>  fAILitail  AUfVB. 

Af  if  fl  _  . 

11  U141  P^i.LiifT 

2‘4V 

(  TO  f  *S5  S*n(  1,0*0 

*  / 

S  tPh0t,2« 

_ _ 

_™0inEXSU'«,  Ef484lin»jrj}StUI 

...  . 

..  . . . 

.  S’P.NOi.i? 

JM* 

t'aylvAwU.CE  IMJAIIU 

IMJel  1  t  1 

StPSOUJ.l 

.Jt*.  .. 

— taulvALtyst.  UiJiin. 

in  jau )  i 

.  *  .  ...  w„. 

StpP'Juil 

J2* 

OAT*  IK|JStll.t»!.Ul 

2’UFtAttO  F*1LU»U  PHui). 

IF  ll  s 

tFSTAS'*-'-  '■'  }! 

J _ LU.AL.Lii A*^  1  liti  i.£LAO  */  ..  _ 

r 

EOvI»*lEvCE  ITIJUU, 

Stlllllll,  M  I2SM  l.lFtjStUI 

S'psaw*' 

JS*.,., 

„CUtt_„  - - 

__  . 

lllt  . .  . 

5t  t  l***Al*^> 

3A* 

C 

S  fi»‘As|VA<l4 

.  E„, 

- r_^ - —  »♦»  ¥  A"  1  *f  Lt  (**l*0T 

£1*1*  *•* 

.  .  5s**»8w2 

C 

6»T*  to  PV  C«A*»,»fl>  TaO*  Tt*K  AA$IC  VALOtS  llstio 

jtPsi.iJ, 

,3»V, 

£ 

...  •*  .4^^  ]h  .ThC.  kOiit 

. S' PS 

>*  oej  i»t*v.  faxM. 

%  iPAllllb,, 

^  i  *,IV 

-t. 

_ ...... _ _ . .  . . .  . 

<  .  .. 

{ 

iO*5S 

V  :  *•*-  V  a  «  i 

»2» 

i 

.  ^..  .  .  •  •  Ml 

(21  li)  . 

,1*1 

%'****'* 

V*  * 

C&*»»(0'fc  F0*T*F  OAi  Ltt, 

«•>  <  0*  » 

.  ^ 

. ?» 'rnoySs 

y  <ni 

I T  >  l*t 

1 

t  i4?7 

%  t  »‘»V,4  '»4 

..  .  £*.*<*.» 

L¥*«i*  .  L^>FR*t 

k'f^4t  -.lv; 

■■  > 

€. 

«*t»CVBta  ST«Us»T« 

STat-sV**, 

*  -  *  . 

-  .  _  Mil  . 

^  i  *  t  MJf 

4  l<> 

4  |C»I 

Sii» 

2  A 

■t'"*  «.  t 

«*t  i 

»*•*%  . 

^  5  **  a  ^ ,  4.  , 

_JL. 

S  t.»»-  > 

« 

i  i  %  s 

1 1  »  1  U*T 

nn 

» iu » 

62V 

■  i  .  -  • 

>5vfto. 

StAKt)* 

i!»*  »».» 

\  4  r  e  •* 

.*'*♦  .. 

C 

.._"..  mt 

.  « j « j  > }  >  1 

>•.!*( 

Sr  4<*t  -  Jl, 

rvAw*  »  * 

2e»»>»  . 

rV4<M5  * 

VJ-»-  -»>* 

.4  / A  f 

18 

V,-*'."..--'’.?  1 

vs* 

C 

„ . ., .... ,  „N,  ^ 

tj-fi  «i*» 

1  ft 

t  it?  i 

-.***  . 

__  . 

. . 

>¥l  ,  .  *><VJ 

...  *r.s  , 

.If yi.  t  .... 

S*P*v4»S» 

t 

T.&Sf  ««Sl2L»> 

S  4»»»-v*te 

i..r  .  ....  .  .  .  • . <  >  M  . 

,  I2»‘l  <041 

$  *  ^  ^  4 

*#* 

" 

>  -4A ! 

j5a,vl.-*2 

JA* 

. .  , . iSH  t‘s»  . . 

V-y.  .4 

* 

s  1 P-  •»  ■» 

ti* 

i... 

-  ..  ?.*?.-?  5f.Mvs. _ 

s* 

♦  %* 

c 

i  i  HVJ  t 

»  »  *Pt-Tt 

Sts,;,,  »f 

IM 

4 

S  1  P“4>  *.» 

WS  A*.  Ut«»is*J|Hu  S  J.A  J 

.i  vt*,x}t* 

stp-v  ,  »» 

107 


AMIN  PROGRAM  (Continued) 


Lines  70  through  115  initialize  the  data.  The  standard  values 
are  those  of  Table  A3-I. 

Lines  116  through  132  read  the  case  number  for  the  first  case  of 
the  run  and  set  the  output  device  coae,  NP.  The  "99  Continue" 
statement  is  the  return  point  for  recycling  and  is  followed  by  the 
read  statement  for  the  case  caption  and  the  addition  of  1  to  the 
previous  case  number. 

Line  133  calls  DECRD  to  input  the  case  data,  which  is  confined  to 
any  changes  from  the  previous  case.  When  the  first  case  data  is 
input,  it  consists  of  changes  to  the  standard  data  (but  must  contain 
one  entry). 

Lines  134  through  138  set  buffers  for  the  plot  routines.  Lines  139 
through  142  set  the  output  format  control,  IP,  according  to  the 
sign  of  the  interval,  OX.  If  DX  is  negative  IP  —  0  and  only  the 
summary  items  are  output,  the  tobies  being  omitted.  If  DX  is 
positive,  IP  '  1  end  the  tables  are  included  in  the  output  (see 
alto  line  207). 


WVi  not  uting  FORTRAN  V,  cord  133  may  be  changed  to 
CALL  DECRD  (UN'rLD) 

end  cord*  518  end  552  in  DECRD  should  then  be  changed  os  described  in 
paro.  (b), 
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a)  MAIN  PROGRAM  (Continued) 


Lines  144  through  169  only  apply  if  the  load  spectrum  ordinates 
are  input.  The  input  values,  TXI,  are  transferred  into  the  load 
spectrum  array  PXL  until  a  value  less  than  0.1  E- 19  is  encountered; 
the  rest  of  the  PXL  array  is  zeroed.  Values  greater  than  unity  are 
set  to  unity  (a  probability  of  unity  represents  certainty  and  greater 
values  have  no  meaning). 

Lines  170-171  ensure  that  all  values  of  X  below  XMIN  are  associated 
with  a  load  spectrum  probability  of  unity,  so  that  XMIN  can  be  set 
at  the  highest  X-value  with  this  probability,  and  the  input  data 
reduced  in  volume. 

The  case  number  and  caption  are  written  (lines  172  through  180), 
followed  by  a  print-out  of  the  data  (as  set  for  the  case),  provided 
that  IP  is  not  zero. 

Lines  200  through  205  form  the  factored  design  load  for  the  case, 
and  if  this  is  less  than  the  maximum  value  previously  encountered  in 
the  run,  retains  the  previous  value. 
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a)  MAIN  PROGRAM  (Continued) 


Line  206  sets  the  maximum  number  of  intervals  ollowed  for  the 
case  (if  different  from  the  standard  value).  If  the  input  value, 
RNB,  is  negative,  the  output  control  is  reset  at  line  207  to  IP  3  -1 
and  represents  a  command  to  print  every  line.  If  RN8  is  posi.ive, 

IP  31 +1  and  the  output  table  is  truncated  as  described  in  PROBF. 

Lines  203  through  232  form  the  number  of  calculation  intervals.  A 
range  from  XMIN  to  twice  D3NLD  is  assumed,  with  a  band  edge 
coinciding  with  UNFLD.  The  resulting  number  of  bands  is  com¬ 
pared  with  the  permitted  number,  NB,  and  if  too  large.  Is  curtailed 
at  the  upper  end.  The  highest  value  of  X  Is  compared  with  the 
highest  test  load  to  be  used  and  the  ouumed  range  increased  if 
necessary,  up  to  the  limit  implied  by  N8, 

Lines  233  through  239  initialize  the  mean  strengsh  arrays,  PSM 
and  PSM2,  and  set  the  X  values,  ensuring  that  a  zero  value  for 
X  is  not  used. 

If  KML  is  !  or  2,  the  load  spectrum  Is  formed  ham  the  input 
properties  of  a  double-family  description  of  the  probability  that 
X  is  the  maximum  load  encountered.  Appendix  II  describes  this 
process,  which  covers  imps  242  through  273. 
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MAIN  PROGRAM  (Continued) 


The  s:  atements  of  lines  276-285  represent  the  formation  of  the 
basic  strength  distribution  as  defined  by  elements  15-20  of  the 
input  data.  If  the  fabrication  variation  is  to  be  superimposed 
(FVARA  not  "zero"),  then  STR  :s  used  for  this  purpose  as  in 
lines  286-292.  The  basic  strength  distribution,  PXS,  is  copied 
into  PSM  and  PSM2  which  ore  then  m^^fied  within  STR.  The 
coefficient  of  variation  of  the  resulting  distribution,  VARS,  is 
used  in  line  293  to  define  tho  intended  mean  strength  by  matching 
the  design  toad  to  a  strength  which  is  SAIL  standard  deviations 
below  the  mean.  The  intended  strength  distribution  properties 
are  then  printed  ot  lines  295-298. 

Data  is  set  for  the  plot  routine  in  iines  299-309.  If  the  'no 
error"  probabilities  of  failure  ore  to  be  printed  (KS  ~  1),  the 
heading  is  written  at  lines  315-316.  The  intended  distribution 
of  mean  strengths,  with  no  error,  contains  unity  for  the  band 
containing  the  intended  meon,  but  is  sere  elsewhere,  os  given 
by  lines  322-331 . 

The  intended  strength  distribution  is  formed  by  STR,  followed  by 
the  use  of  PROSf  to  form  and  write  the  failure  probabilities  end 
reliability.  If  IPiOT  ■“  1,  tf  values  ore  then  plotted  ot  line  338. 
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MAIN  PROGRAM  (Continued) 


Line  339  forms  the  error  function  option,  KE.  The  heading  is 
written  (line  347)  and  stored  for  plotting.  At  lines  355-356, 
the  input  data  PFI  and  PPU2  are  temporarily  held  in  dummy 
storage  so  that  their  meaning  can  be  changed  for  the  double¬ 
family  error  function  (KE  =  4). 

DISCR  is  called  at  line  361;  this  modifies  the  distribution  of 
probable  mean  strength,  PSM,  by  one  of  the  four  error  function 
routines.  If  the  fourth  is  used,  PFI  and  PPU2  are  reset  at  lines 
362-363.  SIR  is  called  at  line  368  to  use  the  revised  PSM  array 
for  the  formation  of  a  new  PXS  an  y,  which  is  then  employed  in 
the  re-evaluation  of  the  failure  probabilities  and  reliability,  using 
PROBF  (line  372).  If  IPLOT  =  1,  the  necessary  data  is  transferred 
to  the  plot  routine  buffers,  at  lines  374-386. 

Line  387  sets  the  test  option,  KT,  and  line  394  sets  the  number  of 
tests,  NT,  When  KT  -  1,  this  implies  NT  tests  to  each  non-zero 
test  factor;  when  KT  =  2  or  3,  this  implies  NT  tests,  one  to  each 
of  the  NT  test  factors  input. 

Lines  395-396  duplicate  the  cumulative  probability  of  a  given 
strength  in  PRS2,  to  enable  the  original  values  to  be  retained 
in  PRS  for  later  use. 
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a)  MAIN  PROGRAM  (Continued) 


The  operations  on  this  page  incorporate  the  effects  of  NT  tests 
surviving  each  of  the  input  test  load  levels,  defined  by  a  test 
factor  applied  to  the  un factored  load. 

The  headings  are  formed  (and  stored  for  the  plot  routine)  in 
lines  400  through  410.  The  counter,  MT,  is  initialized  as  1 
and  the  first  test  factor  tested  for  a  "non-zero"  value.  Headings 
are  written  (and  stored)  in  lines  417  through  437,  and  TEST  is 
called  to  form  the  probability  of  surviving  all  subsequent  tests. 
6AYES  is  then  used  to  update  the  mean  strength  distribution,  PSM2, 
this  then  being  used  to  update  the  individual  strength  distribution, 
PXS,  by  means  of  STR.  PROBF  is  used  to  form  the  failure  distri¬ 
bution  and  reliability  and  the  values  are  stored  in  the  plot  routine 
buffers  (lines  445  through  457). 

This  process  is  repeated  for  the  remaining  tests  to  the  first  test  level 
(the  loop  from  line  423  to  line  459).  MT  is  increased  to  2  ond  ths 
whole  process  repeated  by  a  return  to  line  412.  When  a  "zero" 
test  faetor  is  encountered,  the  case  is  ended  by  returning  to  line 
129. 
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MAIN  PROGRAM  (Concluded) 


The  final  operations  are  similar  to  those  on  the  previous  page,  but 
are  for  the  remaining  two  test  options.  In  both  cases,  NT  tests  are 
made,  one  to  each  of  the  NT  test  factors.  When  KT  =  2,  the  tests 
result  in  failure  at  a  load  assumed  to  lie  In  the  band  containing  the 
value;  for  KT  =  3,  the  tost  survives  the  specified  load. 

The  appropriate  heading  is  written  by  line  463  or  468  and  is  stored 
for  plotting  in  lines  47V  through  477.  The  pre-test  distribution  of 
mean  strength  is  copied  into  PSM2  which  is  then  updated  by  BAYES 
(line  497),  and  used  to  update  the  individual  strength  distribution 
(STR  is  called  at  line  499);  these  new  values  of  PXS  ore  then  used 
in  PR08F  to  re-evaluate  the  failure  distribution  and  reliability,  which 
are  stored  in  the  plot  buffers  (lines  503  through  514), 

This  loop  (lines  480  through  515}  is  repeated  for  each  of  the  tests, 
after  which  the  program  returns  control  to  line  129  for  the  next  case. 
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b)  DECRD 

This  decimal  read  routine  is  similar*  in  effect  to  the  one  used 
(but  (  at  listed)  in  reference  1.  The  description  in  reference  l 
remains  applicable  and  is  reprodi  ?ed  here. 

In  the  Decimal  Read  data  input  method,  each  card  is  divided 
into  six  fields,  each  containing  12  columns.  The  first  field  is 
reserved  for  the  index  which  is  the  Data  array  location  of  the 
data  in  the  second  field  on  the  card,  so  that  five  fields  are 
available  for  data.  However,  it  is  not  necessary  to  supply  a 
number  in  each  field;  if  a  field  is  blank,  the  program  will  retain 
the  variable  unchanged  from  the  value  already  stored  in  the  DATA 
array.  The  remaining  four  fields  on  each  card  represent  the  loca¬ 
tion  of  variables  which  ore  in  numerical  sequence  after  the  first 
location. 

c)  EVP 

This  routine  forms  the  distribution  properties  of  the  load  or  strength, 
using  Gumbel  equations  for  double-families.  For  the  load  spectrum, 
maximum  extremes  (K  <=  +1)  ar*  used,  the  tall  extending  towards 
higher  loads;  for  the  strength  spectrum,  minimum  extremes  (K  *  -1) 
are  used,  the  toil  extending  towards  lower  strengths. 

Lines  1  through  10  allot  storage,  etc.,  and  are  followed  by  definition 
of  the  intercepts  ond  slopes  of  the  characteristic  straight  lines  of  the 
Gumbel  plots.  The  overall  mean  of  the  double  family  is  formed  at 
line  30,  and  the  summation  terms  ore  initlolized  in  lines  31-34. 


*WH«n  not  using  FORTRAN  V,  cards  518  ond  $52  should  be  changed  to 

SUBROUTINE  DECRD  (RK 1) 

550  5TOP 
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c)  EVP  (Concluded) 

The  basic  loop  comprises  lines  35  through  85.  For  each  band 
of  the  variable,  X,  the  cumulative  probabilities  are  formed  at 
the  band  edges  for  both  of  the  families.  The  differences  give 
the  population  in  the  band  as  PHT  (line  68),  double  precision 
being  employed  to  improve  accuracy.  First  and  second  moments 
are  found  assuming  the  band  contents  to  lie  at  the  band  center. 

For  the  load  calculations,  the  resulting  values  are  integrated  to 
give  PXL,  the  probability  of  exceeding  X;  the  integration  is 
actually  formed  by  subtracting  successive  increments  from  an 
initial  value  of  unity. 

For  the  strength  calculations,  the  increments  are  stored  in  PXS 
and  the  cumulative  probabilities  stored  in  PRS. 

Finally,  the  mean  (BART),  standard  deviation  (ST)  and  coefficient 
of  variation  (VAR)  are  formed,  followed  by  the  return  statement  at 
line  89, 

Lines  90  through  96  are  provided  to  remove  a  possible  anomaly  when 
family  B  is  to  be  subtracted.  If  the  cumulative  probability  reduces 
in  passing  from  one  band  to  the  next,  a  negative  population  is  im¬ 
plied  for  *hat  band,  which  is  physically  absurd.  The  presence  of 
this  irrational  value  is  detected  by  a  near  zero  increase  in  the  trans¬ 
formed  probability,  YT  (line  59),  When  such  a  condition  occurs, 
BETAB  is  raised  by  five  per  cent  and  a  new  attempt  mode;  this  is 
repeated  until  valid  results  are  obtained. 
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d)  STR 


This  routine  evaluates  the  resultant  strength  distribution,  PXS, 
and  its  cumulative  probability,  PRS,  given  a  distribution  of 
mean  strengths,  PSM2,  and  a  basic  definition  of  the  shape  of 
a  distribution  with  a  given  mean. 

Lines  I  through  7  define  the  basic  storage,  etc.  and  the 
initial  values  of  the  working  arrays  PXS2  and  PRS2. 

The  loop  of  lines  8  through  30  takes  each  mean  strength  level 
in  turn  and  forms  the  distribution  of  that  contribution  to  the 
total,  using  EVP.  The  inner  loop  of  lines  21  through  29  sums 
the  contributions  of  each  sub-distribution. 

The  second  phase,  lines  32  through  57,  permits  the  superimposition 
of  a  second  variation  such  as  that  due  to  fabrication,  and  reforms 
the  values  of  PX52  and  PRS2. 

The  remaining  lines  sum  th#  total  and  the  first  ond  second  moments, 
farm  the  overall  mean  ($ART),  sfcndord  deviation  {ST)  and  coefficient 
of  variation  (VART)  and  also  copy  the  working  arrays  into  the  common 
arrays,  PXS  and  PRS.  before  returning. 
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PROBF 


This  routine  calculates  the  incremental  probability  of  failure  (DbLPF) 
for  each  strength  level  and  integrates  to  give  the  cumulative  failure 
risk,  PF,  The  final  volue  of  PF  is  subtracted  from  unity  to  fctm  the 
reliability.  The  routine  is  also  used  to  print  oil  relevant  output  values. 

Dimensions,  etc.  are  specified  in  lines  1  through  20.  If  the  output 
code,  IP,  is  not  zero,  the  heading  is  printed  (line  22).  If  IP  is 
negative,  the  start  control,  JST,  is  set  to  one,  otherwise  it  is  set 
to  the  lost  interval  hcving  a  load  probability  of  unity  {lines  26-32). 

For  the  first  interval  (which  will  always  have  unity  for  load  prob¬ 
ability)  the  cumulative  strength  probability  is  used  to  form  the  total 
probability  of  failure  at  this  load;  the  values  ore  stored  for  plotting 
and  ore  wri'  sn  at  line  45. 

The  remaining  intervals  are  then  treated  in  turn  in  the  loop  (linos 
4p  through  74).  When  IP  is  ♦'I,  the  insignificant  lines  ore  not 
printed;  th^se  ere  chosen  as  those  for  which  the  incremental  failure 
probability  is  less  than  KT^  and  for  which  the  probability  of  o 
lesser  stren^'h  exceeds  0.999995.  The  lost  line  is  always  printed. 

The  resultant  reliability  is  set  by  line  73,  load  levels  below  yn- 
factared  load  being  jumped  . 
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_ QlflEtL5J.Qii-l.-8y  EiSPi _ 

EQUIVALENCE  { EauK ( | )  iLBUF lit) 

_ DOUBLE  PRFOlSIQN  0  HE _ _ _ _ 

f  •  •  • «  F 

_£ _ jq«.  aarifi-tt-puas _ _ ___ _ ___ 

0 1  HENS  I  )5l  BUFR  (  12)  ,H0SC(1»  )  .KI24I3)  ,F|32I  1 3 1  i 

_ ! _ K22  l»»'(|pLlN|PIBltriLIN2Pl8l  .HL1N3PU)  iHLINflPItl _ 

iF'ijp)  If, 21. 20 

20  iiLUXIMiLPi _ 

10  FORM  A 1  1/Jl  X  ,  '  t  *  I*x.*x*  .BX^PXL'  .«X,*PXS'  .BX.'PRS'  ,  7X  ,  *  OELPF  •  ■ 

_ L8-L  i-'J'-E-Lt-I  A  i.  jjPSflJLl _ 

21  ONE" 1 • 

_ ifi  IPJL  25 ,24x26 _ _ _ _ _ _ _ 

25  JST"1 

_ TXL"1  . _ _ _ _ _ 

<50  TO  25 

_ 26  JST«NSt  _ _ 

TXL«PXLf  .1ST  I 

_ l<L.IS.T«J5Tti  . . . . . . . .  . . . . . 

OELPF"TXL"PRSljSTl 

_ PF«OElPF _ _ _ 

IF  <  IPLOT  .E<5.  0)  <50  TO  RUOO 

_ XPI  TTIJUXUSJ) _ _ _ __ _ _ _ 

APLTT(2l«Xl JST) 

_ »PLT Bl 1 I«X1 JST I _ _ 

YPLTTl 1 ) "OELPF 

_ YPLTdl I 1«PF _ 

IPLT.2 

_ _ k!-L _ 

6000  CONTI  NOE 

_ IKI  IPlJ?Ji2Jt22  _  _ _  _ 

21  AMI Tt <NP, 11  I  JST,* ( JJT) .TXt.PXSt JST ) ,PrS( JST ) ,uE»>F ,PF ,PSm2 (JST I 

22  00  1  1  " 1  ST .NX _ 

I K  < I -NS  T I  22,28,28 

_ 22_0ELPF*',|tSl  1.1 _ _ _ _ _ _ _ 

T<L"t»0 

_ ...GO  fO_29  . .  _ _ 

2»*Ul>>*»i(  l) 

_ Oi:LPF"TXL»PX5t  I  ) _ 

25  PF"PFiQELPF 

_ if  cploj  .lo.  oi  ao_.ro  idoi _ __ _ _  _ 

L"L  ♦  1  .  . .  . 

_  YPLTtJ (l I "PF 

ypltt< iplti"oelpf . 

_ |PLT«  1PLTH _ 

«Pl  T  T  I  I PLT 1 *X  t I  ) 

_ YP,L,.T  T  (  I  “l  T  I  “OtL-PF _ _ _ _ , _ _ _ 

I  PL  <  •  I  PL  T  4  1  . " . . 

_ , . XPlTI  1JPL.T1".X.U  1-- _ _ _ _ 

APLrO<Ll"*( I  I 

6'lQI  CUNT  i»UE _ 

.  i  r  ( i  p  i  r  ,  j ,  b 

_ S  IF<UELPF-«te-lO)  6,6,1  _ . _  _ 

6  IF  IUM-PRSI  I  I*.1F»S)  7,7,V  "" 

1M!'-N»I  2.1.1 

1  JHlTEtNP.il  J  l.Xlll.TAL.PXSlli,PHSlll»OtL  F.PF.PsniUl . 

11  FORMAT  <  1 X  , !  2  «f  1 '  -  A  **C  1 1  •  t>  2 

.  3  u.iw  -  - - ...  — .  . 

21  ARL"0NC»PF 


STPRU721 
_S_TPR0725 
1  5  T  PH  0?26 
5  Tp  ho  ’ 
STPHU728 
STpRU729 
ST  PR0720 
STPROTo 
STPhJ/22 
„Sl£ilOZ28. 

ST  PR0721 
_S.TPitO.7Q5 
STPK07J6 
STPR0727 
STPH0728 
5TPR07 J9 
~STPR07l0 
STPR0711 
STPR0712 
STPKO/U 
STPR0711 
STPRQ71S 
STPR0716 
STPRQ717 
STPR0718 
5TPR0716 
ST  PK07SQ 
ST  PRO?  5  I 
ST  PH0752 
STPRQ752 
S1PNU7S1 
5 1 P  ^ 1755 
_STPPu75i 
STPK0757 
SIPRU750 
STPN0759 
"  S T  PR07  6  0 
_5  T  P  R  0  7  6  ) 
S  T  P  R  O  7  6*2 
5TPNU/62 
ST  PRO  7p1 
_S  T  P  KO  7.6  5 
ST  >'«  07  06 
5TPR076  7 
'  S  T  P  R  ( J  7  6  M 
5TPR0769 
STPR0770 
STPR077I 
ST  PR  j7  7  2 
STPm0772 
'  ST  PRO  /  79 
ST PK  0  7  75 
STPR0776 
S' PK07  7  7 
S' PRU7/R 
STPR0779 
STPmO/RO 
S 1  P  N  -j  7  n  I 
Si  P R 0 7 m 2 
_  S  T  P  R  0  7  d  2 
STpRu7ui 
STPWO/tiS 
”STPRo7Bn 
5  T  P  R  ■?  7  u  7 
S  T  P  R  0  7  6  ti 
STPR0769 

’siprj'/vu 
_AJ  PJiQZta. 
STPRU7V2 
.5TPR07V2 
STPRO/99 
SlPhJTVS 
STPRJ79* 
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PROBF  (Concluded) 


The  Final  Sines  print  the  total  failure  probability  and  reliability, 
storing  the  values  for  plotting  before  returning. 

BLKDAT 

This  subroutine  stores  the  titles  used  in  the  plot  routine. 


.71* 
7S* 
7  4* 
7/* 
73* 
7V* 

41* 

b<* 

Si* 

gM* 

«4* 

.84* 

s7« 

as* 


-S-LPi-L/ii. 


MITE  lrfP.l2l  AKl.PF 
12  I  u'KMAT  1 SX » 

2  ‘asymptotic  reliability 

.3*  PRUB.  £13.7/1 

If  UPLOT  .£«.  Cl  f.O  to  SCO* 
ENCOflb.l.lEltUF  .jUlilj-LtK L _ 


INDEX  IS  *.F  m.7,SX,'F*lLuHL»  , 


HL  II.aPlSI'F.BC*- 1  } ) 

.  Hut  HIP  141 •EttCl t  2 1 
,u  |  NHP  t  7  1  *£  HUl  (31 
HLihtPl  A) «t3wF  t  *f  J 
S.TO?  FORMAT  1  f  1 0  *  7 1 
SQQ2  CUN 7 1  :.u£i _ 


RETURN 
EwO 


Si  P*;i7v.i 

stpho/vy 

STPPJ6U0 
StPHJbu J 
St  PKOtllli 
S  t  PR'xn)3 
STpRobcn 

STPPutiCS 
S  IPHOH'Jt 
5tP'<SCw7 
sipKca^a 
St  P  R  J  0  u  ‘1 
StPHJbl J 
StPKObl 1 


1* 

2* 


bloc*  n*T* 


STpKum 1 2 

•>t  >>h,ih  I  -I 
S<  PKOb  1  S' 
StpSJblS 
S  T  p  R  f;  f  |  a 
SlPhJBI J 
Stl'hjnis 
5  T  P  f»  C  b  1  V 


3* 

a» 

s* 

4* 

;• 

_j_*_ 

v~*~ 

10* 

it* 

12* 

13* 

n* 


C  OLliciC  0  A  T  *  SUPFROsRAM  to  contain  titlcs 

{  FOR  CACCOMP  PLOTS 

COMMON  /COmPlT/  XPLTTIS02I  •  VPuTT  ( 102  )  . T PL  To  12021  ,aPlTa<2i)2) 

_  1  ,oUFR.H0<*C.FU4.f  132.F22.Hl  t  N.  1  P  »HL  1 N2P  . 

2  HUnJP.HLlt.HP 

QIHUiSlO'l  SUFR  (  I  2  1  .  H  )SC  (  IT  )  .M24t3  1  .Fi32{  |3)  , _ 


l  f  22tb7l  ,,(U!«lr'(8l  ,rtLIN2Pt#l  .hLllliPtSI  *PU*iSP18) 

1)1  HENS  I  ON  E  22 1  V  1  . . 

EQUIVALENCE  1F221 1 1 .E221 1 1 1 
OATA  l»06Ct I  I . t«t ,3 )/*CASL  SO.  _V 

’  OATA  1 124  /’TEST  SERIES  •/ 

Data  f  132  /’TEST  NO. _ TEST  ftClPN _ 


■  TS.s.T.  LOAD 


STPAUPaJ 
STPhOb/l 
STPf 00?2 
StPKOb/J 
STPHjS/A 
STPRU.I2S 


IS* 

14* 

1 7  * 
i«* 
IT*' 
20* 


•/ 

/«USFACTO«EO_LOAO  ■ 

/ •UnOEVs T RENT}  f H  RISE  '  *  . 

/•RELIABILITY  index  «  _ 

/•ASYMPTOTIC  *EL.  index  * 

E  23  /'PKOB.  or  SURV I Y 1  NO  NEAT  TgSTS  TEST  lOAP 


l 

OATA  RLlNlP 
"  6aTA  ht I N2P 
OATA  HL1N3P 

oata'hlVnsp 

OATA 

T7 - 

End 


•/ 

•/ 

”""•/ 

PHOB. 


StPKUhib 
STPb'JPiT 
SI  PHI)l*2b 
ST  PKOBW 
STP,Hu'o30 
•STPKObO 1 


21* 

22* 


ST  PA'J«J2 
STPbOU.*3 
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9,  SUBPLT 

This  is  used  in  conjunction  with  ti  .  system  routines  PLOT, 
PLOTS,  SCALE,  CAXIS,  SYMBOL  and  LINE,  to  plot  the 
appropriate  output.  Values  are  written  cn  magnetic  tape 
for  eventual  preparation  of  hard-copy  output. 


SOOMOOUnE  SUOPlT  lf**LT,MXI  STPKOoiS 

COMON  /COPPLT/  HPlTT  C  H02  I  .tPLTT  MCi2i  i  VHJU  Tt}t.2Ui  )  ,XPLTBI2U2)  ..  STphMtJS 

1  ,00FRtH0<iC.F12*iF  IJ2.F2*.hLI«|P.MLlN2P,  StPi-o.‘J* 

2  KLlh.JPtRl  IQHP.FftUFlSJ) _ STFtuui.I/ 


s 

6 

/ 

a 

» 

_Li 
1 1 
12 
i  a 

_  1H 

is 

_Lt 

1/ 

I* 

'  1» 
_iU. 
2 1 
22 
2J 
2H 
2S 

2? 

2a 

2» 

JO 

)1 

-ii. 

JJ 

Jh 

J4 

J* 

J/ 

JL1 
Jy 
HO 
HI 
Hi 
H  J 

a  h 

Hb 

'1* 

'hf 

na 

"hv 

&u 

si 

V2 

Si 

%H 

S4 

'■i  6 

s?' 

SO 

5» 

AU 

a  i 

Ai 

~Vj 

*H 

a  l 

*7 

*4 


0 !  HENS  1  On  bUM  I  I  3  »  .HltfiC  1 1 9  )  ,  F  1 26 1  J  .* ,  fl  J2 1 1  J )  t 
t.  F22IA»!  |Ht,lMlK(8)  .HUINiPtH)  ,HLIM3P'«.I  .HI.INHPIB) 

0 (PENSION  bUFFEKI 102HI «  ICESTIST 

DOI  I0ENT  /  •  SCj7»t700  07-IN  bRTObf  .  ‘/ibLANK 

Data  mis  /<»•/,  YAxtsa  /•psfj*/.  yaaist  ’/’opif  »"•/ 
if  1 1 Plot  ilti  qi  oo  ro  7  a  a  u _ 


•  / 


sippoaja 

STPnOBJV 
sTp*<obhu 
stpkobh  i 
STphOhhX 
SIPRUbHJ 


f  •••< 

,t _  iNlTHUit  PLOT  .  PACKAISC . . 

MM 

_  C«LL  PLOTS  iBUFFIH.HUItlwa.tOrWT.JO) 

LPtOT*-l 
7r3u  CO.NllM.iE 


S IPKOOhh 
SlPRUbHS 

S I PhObH A 
Siphons? 
STpHuoNrt 
SlPhliat? 


.  <•*.* 
C 


call  Plot  io.o,o.o,»2> 


CiE  1  H|N»  A  ANO  DEL 1 A  A  •  SAflE  FOR  TfP  AND  OOT  IOK  P|01b 
Call  SUiX-U£iJUL>JlJLlUJULaJ _ 


STPhOtbU 
STPWIjnb  l 
‘S1ph\.aS2 
JiLLPMAttii 


IlCTOP*N«*2 
_ APLU I  JRCTOP*H«XPLTfllN*»U _ 

0LPLT»APLTBtNA*2) 

_ XPLI1.I  LMCTUP*2)«XPLTa.lN.AliJ._.. 

XHNPL*APLT8(NX«1 I 


STphutibR 

STPIOhbS 

STPh,H>!>* 

_S.tpkobs; 

STpMuuva 

SIPROnSY 


MIN*  r  A  NO  OClTa  y  FOR  UOTTOH  ANO  TOP  PLOTS 

_.call  scale.  I.ypli.R»1».5jN*».*.I _ _ _ 

Yl1l'(Pt.»»PLTlltNX*l  I 

Ore»YPLTR( N«*2  I  . . 

CALL  SCALE  I »PLTT i J.Si«X»2, I  I 

YIIPlPTnYPLTTI  InCT0P«I  I _ 


STPNOBou 

STPhUUAl 

STPWOltnX 

STMhOHftJ 

~STPRO*>«h 

51 PRJnnS 


OYY«YPlTTI  i Nc 7  0P*a 1 
.CALL  L  A. «  i  S  1 0 . 0 . 0 . 9  ,  X  a  <  I S  »  - 1  » (|  •  0  1 1 . 0  (  A  «J.NPfc  i  OX  P  L 1  Ol  •  I  Uj  -  J 
CAIL  f  AA|S<OtO",l«btYAAlSR,R*jAS»*U*UtYt«|  AP0*0»ti.LUlll«2> 

.  C  ALL.  S  Y  H56l  «  J.Sil. 1.0*10, HUMfA«,»*li(U.. 
if  i  |p|.T  ,cw«  iii  i»o  ro  20 
call  SYHSOL  I  I  «  0  i  1  tf  1 1  t  0 «  1 0  if  i  37  1 0  ill  i  »?>  t._  . 


SIPNOn** 
STPKyun? 
'STPhiau.li 
,STPhOh*R 
SlP»Uu7U 
5ipRi;a7 1 


Call  syhaol 
..  CALL  SYMaOL 
OVNHJ. J 
_  Y  *a  t  v 

’  n«u 

_ ao  IS  l«l.lPLt 


(StS.Y.4,U*lilir22il)t0iJ0) 
IH./|V,2.0,IJ,F22I*I  AOh.fiilt}..., 


CALL  SYHAOL 

t  j*  U»V  _  , 

V* Y  *0  >  R 

IS^COUTmoE _ 

20  COPT (HOC 


(H,7,Y,U.10.F22UJI  .D.U.JOl 


-CJJU 


C  *  •  «  • 

<  ..  . 


ia*c  hoy  ton  plot 

CALLLjNt  IXPlTHiYPlY.I.NA.I  |0»2*I  . 

put  axis  on  top  plot  amp  oraa  it 

CALL  CAAISIO*OiS.U.AAAIS.“l  i  R  *  0  t  0  •  U  .  AM  I  hPL  »  OaPl  T  .  UlOi*U 
call  LA<IS<U<IJ  »S.r.  YA<IST(1mJ»S»7U.O.YaI»»T(dYT,U.1J.2I 


ST  PROS  72 
SlPROS/J 
S  T  PRO*  7  H 
SlPhpbfS 
STPhbb/o 

5JL?  **  s  ajif. 

■$1  PR OH  7 8 
SlphgtsfY 
5  TPnUnwO 
STPRvipal 

S»PMii«ti2 
_S  T  PhOnn  J_ 
“sTphobn™ 
SlPhoPhS 
si  PH..8U* 

.  STphdH«; 
si  Pht'hpa 
SlPhQnoY 
Sfpii'ytt'vO 
STPhoivI 
STPhteHA' 
SJpp.iavJ 

STPNuoVH 

s  t  p  “  a  *  r  s 

“sTpViTsh  *" 
STPn,iu  »7 
sTphviaao 
jTp»ufc»» 
St  PhbVnv) 
STpnOVuJ 


fiat 
C  _ 


f  •  •  »  • 
C 


PuT  6«  I'ACiE  mlaO|..(*S 

CALL  SYHROL  H  .0. 10*5*0*  Ib.HniC  ,Ut{l,nai 
CALL  SYOSOl  I  I .Ut IOi JtQ» lOihOLC >»• iO*y*PAI 

call  syaaiL  (j»t>>v»R,o«il|»(««FAiit»o«o*J4* 

CALL  SY.HHOl  I  J.  5  I  V  «  t,  It  I  ftitlUFH  <7  >  iQidt  >A  > 


RfORISlN  FOR  TjP  PlOT 
CALL  PlOMU«UiViUi*J) 

Call  Lt.N.e  A  APL  t  T  f.2 1 1 YPLT  T 1 2  I  •  I RCT  OP^*  I  *  1 1  S»1V  1 

RtTORh 

EhO 


DISCR 


This  routine  uses  the  appropriate  error  functions  to  modify  the 
mean  strength  distribution  array,  PSM.  It  contains  four  alter¬ 
natives,  described  in  Appendix  il. 

Lines  1  through  7  set  the  storage,  etc .  anU  the  program  then 
splits  four  ways. 

If  K£  is  1 ,  the  Bouton/Jablecki  function  is  formed  and  used  to 
define  the  probable  PSM  values,  as  implied  by  lines  9  through 
23. 

When  KE  is  2,  the  Freudenthal  function  is  used  in  the  same  way 
(lines  24  through  37). 

For  a  KE  of  3,  a  Gumbel  function  is  used  (lines  38  through  51)* 
The  fourth  option  differs  in  using  a  double-family  distribution 
whose  means  and  coefficients  of  variation  are  input.  Lines  52 
through  72  are  similar  to  the  corresponding  ports  of  EVP  tmd  also 
form  a  modified  PSM  array. 


jmmvUKS.-atsc1*  LS£jjmaj»uiej 


frlPK  ;<i  i 


2*  COMMON  (/*/  f*  HU  I  200  t • A 1 2QQ  >  ,t.A  .MS1200»  ,MS  1 2L0 »  (PS" I  2[)0  1  ,PSM2t20G)STP*0Vul 

_J* _  COMHQt./pAl>/UNKL0.rt.01.4KH.HHML.Ln*HA.L  ■' *H  A  ■  L SU«d  .  I B A»U  .  I V  AH B  .  ...  SlPMOVClt 

*•  fiAl.L.«Sl8SNWO.MAS.R«MSiS8*HA.SVAKA»S$Ort«,l.rtAil».SVAR»|FhAHA«FVA«A.  SlPKORUi, 

-A! - -iF.jjUHttAiaARO.FVARB.HitF.pn  .PPlH  .PF2.PFU2.HICT.HHT.UHH _ 5LlEJ1U.RU'-. 


A  * 

<» 

io»' 

JUi 


UCAL  LaARA,LVAHA,LSUM*,l.8AUb,LVA«B,MS 
_.BM.UbLE_J’i1EC  1 5  I  ON.  M I  .<*«2,i-a |  ,Pj|2  ,t*$  I.,«>s2 
00  TO  ll.Otl . I&) ,  <C 


_2_Atio»iv‘jiAU05io(Pf:iiMi.QaioiFFuu*Ava<ai'J?Fl.LiAkou»iaF.Eujti)/ 
I  I  ALOetO(FF'2)«Av(VG)')IFFt>  »  1 
»a«i,/a 


STphOVk/ 
SIphQVIjB 
SlpHiiSu* 
SlKMOl  | .1 
STPntjti  I 
SlPH JV )  2 


12* 

_L4*__ 

n* 

JS* 

'»*• 

JLit, 


«  !•(  ALOOIO  (PF2  l-ALOOLatPA  l  I  I  1  ULQU 1 0  IPPU2  l-*i04  1  0  (  PPU  l  I ) 
J»_O..R  l  •.!,>»* _ _ _ _ _ . _ , _ _ 


PHl*(H**(*t  |  )».SAOAI/OSWLUi**l»| 

_ a..i»i  ...  . 

I  >M2*<NA»|«m  *.S*0»  l/OSULU)  I 
IF  I  P  H 1 «  )  .m  11,13.11 


SlPKI'Vl  i 
s  Jpwtiv  i  •» 

SIpHOUt. 
S  Ip*  JX  i  o 
S  I  P«H,U  l  > 
Si  x nu>  l  n 


ti*  10  P<2«|.0 

_ it  .. 

*q»  iMPsnm-.te-iai  *.s.<* 

21*  FSrtlt-O.O^  . . 

ti*  R  CONtJUJ? 

?>«  -  <■->  to  j 


in* 
F»*  , 
n* 
4 *' 
21* 
JU> 


0  SI  *Ui  OSkol- At  00  l  l.0-*F2l  >-Al.ft4tOI-AwOG«  UO«l*ri  1 1  1/ 

_ _ )  i  tuSfiV-i*  >*u?  >»iuo<(  1 0 1 -'Flu  *  I 

UA>H«4UO0l  t  ♦0*»»  lll**»t»/ajl  t  /PPvkl 

BO  *.  . 

PHl»5APi"<H**llti  l«4S»{lt  1 10S«4D»***l  > 


40. 

.31* 

»i* 

J*. 

»• 

14* 

«0» 

'  I  * 


*  UA*lll  I  1  *.%*0A  1/CSvi.ol  **M  J 

...  ItlKM.OI  IN.It, 1» 

it  s*i*i,a 

J.1*.  *tSM  1 

1»m  i  |  »«0.3 _ 

— aJu  *  I  '■  tvin.  ■  ■  r4'  4  -m  >W—  mm 


ni* 

mi* 

«** 

«»• 

*»* 

,  «?* 

»t* 

4-1* 

;.  41'* 
%2*- 
Si* 

’s*»“ 

/*** 

s** 

42* 

**» 

■J1‘a 

*4* 

Hi* 

„>3*  . 
**« 
*S» 

♦  ** 
*/• 
*.** 
«♦* 
10* 
11* 
ti* 


•  COhMnuC 

Oil  TO  JO  .  „  .  _ 

AI»*V.Oil  iSSUcOSI  I  «3»^f  till 
*4»Atl<jc,i4*si»to&n  ,a«*rJs i  i 

\  U|„U»i I  *»t*j|* t 22 1 *(*>>. lil (*»*»!  > 


OJ  »0  IM*** 

<l**Sm2M-2«2SS,M*M*l*'<M  t  I  *  *4*0* » 20*41.0  >/S 
1 2“t  I  “I »  SS<SS*3*  i  i  tts4vft*S 1 
*«l*St.*lTl  | 

PH«*4iVt »J| 

>S»  Itl  «0*S>  t»'jA|.«l>«2 1  _ _  . 
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h)  01  SCR  (Concluded) 


¥ 


The  final  lines  are  common  to  all  four  options  and  are  a  test 
that  virtually  all  of  the  distribution  has  been  formed  within  the 
permitted  X-range. 

5)  TEST 

This  routine  evaluates  the  chances  of  surviving  each  of  the  subse¬ 
quent  tests  in  a  series.  Lines  1-20  set  the  storage  tocottons  and 
write  the  heading. 

For  NT  tests  to  each  specified  (non-zero)  test  factor,  lines  22 
through  50  locate  the  band  containing  the  test  toad,  XT  and 
extract  the  probability  of  strength  less  than  this.  The  complement 
gives  the  required  chance  which  is  self-multiplied  for  each  test 
in  the  series  (line  33).  Values  are  written  ond  stored  for  the  plot 
routine  before  reti  rung  ot  line  50. 
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TEST  (Concluded) 

Lines  51  through  78  perform  a  similar  function  for  the  alternative 
options  of  one  test  to  each  of  NT  different  test  loads. 

GEV 

This  function  evaluates  F  =  exp  (-exp(-y))  for  the  Gumbel  dis¬ 
tribution.  Double  precision  is  used,  and  four  ranges  ofy  are 
separated.  For  y  <  -4.0,  the  probability  of  a  lesser  value 
is  negligible.  For  values  between  -4.0  and  14.0,  the  basic 
function  is  computed.  When  y  exceeds  30,  the  probability  has 
reached  an  effective  limit  of  unity. 

For  values  between  Hand  30,  direct  computation  breaks  down 
since  exp  (-y)  becomes  insignificant.  Expanding  gives 

-  1  e-y  +  e~fy  +  e  — —  A3- 1 

21  ~  31 

but  since  y  is  large,  all  but  the  first  two  terms  tend  to  zero, 
and  the  function  can  be  replaced  by  l-exp(-y). 
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BAYES 


Subroutine  BAYES  uses  Bayesian  techniques  to  update  the  mean 
strength  distribution,  PSM2,  so  os  to  incorporate  the  test  results. 

Lines  1  through  18  allocate  storage,  etc.  and  initialize  volues. 

The  basic  loop  (lines  19  through  45)  first  forms  the  constants  of 
the  double-family  distribution  of  strength  with  mean  at  x(i). 

For  survival  tests  (KT  3  1),  the  probability  of  a  value  greater 
than  Xy  is  formed  at  lines  25-26  and  32-36.  For  failure  tests 
(KT  =  2),  the  probability  that  the  test  load  is  in  the  interval 
x  +  l/2dx  is  formed  by  lines  28-36  and  38-41 . 

Line  42  forms  the  numerator  of  the  fraction  on  the  right  hand 
side  of  equations  A2-48  or  A2-50,  the  PSXT  value  being  appro¬ 
priate  to  the  type  of  test.  Lines  43  and  44  determine  the  maximum 
value  of  the  numerators  for  ail  bonds,  and  line  45  sums  the  volues 
to  set  the  denominator. 

If  the  denominator  is  “zero'* ,  lines  47-50  set  o  diagnostic  message, 
and,  if  the  summation  consists  essentially  of  on©  term,  lines  51-56 
print  an  appropriate  warning .  If  o  valid  expression  exists,  lines 
57-59  form  the  posterior  distribution  of  mean  strength  os  PSM2. 
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A3. 5  User's  Guide 


a)  Deck  Set-up 

The  first  two  cards  required  are: 

@  RUN  Card  containing  run  identification 

#  LOG  Card  giving  accounting  information 

where  the  symbol  (3*  represents  the  multiple  7-8  punch  in 
column  1 . 

The  FORTRAN  cards  for  the  source  decks  of  the  main  program 
STRP  and  the  ten  subroutines  (DECRD,  EVP,  STR,  PROBF, 

BLKDAT,  SUBPIT,  OISCR,  TEST,  GEV  ond  BAYES)  follow, 
each  being  preceded  by  on  input  and  compile  cord  of  the  form. 

(«  FOR,  IS  .STRP,  .STRP 
After  the  !<nt  routine,  a  cord  is  inserted  bearing 
.«  XQT 

end,  if  plots  ure  to  be  generated,  this  is  preceded  by 

#  ASG,  T  4.,  T,  reel  no. 

b)  Data  Input 

The  format  has  been  described  in  paragraph  A3*4fb),  and  0 
sample  input  set  it  shown  in  Table  XXXV .  The  first  cord  con¬ 
tains  the  cose  number  of  the  first  cose  in  the  run;  this  is  right 
justified  to  Column  S;  if  output  device  10  is  to  be  used,  the 
case  numb«r  is  negative. 

The  next  cord  bears  the  coption  for  the  first  cose  (op  to  72 
columns,  12A6  *nrmot  being  used) .  This  cord  must  be  present. 

The  remaining  cose  cords  bear  input  data  which  differs  from 
the  standard  built-in- values.  Each  cord  carries  up  to  five  values; 
the  location  of  the  first  is  right-adjusted  to  column  12  and  the 
other  locations  ore  implied  as  consecutive.  The  arithmetic  values 
or*  EI2  former;  decimal  values  rocy  He  anywhere  in  the  12- 
column  Held,  but  integers  ond  integer  exponents  mutt  be  right 
adjusted  to  columns  24,  36,  48,  60  or  72.  th  lest  data  cord  for 
the  case  must  hove  o  minus  sign  (-)  in  column  1 . 
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TABLE  XXXV 


A  similar  set  of  data  follows  for  the  second  case  (the  caption, 
followed  by  any  data  changes,  the  last  card  having  a  minus  sign 
in  column  1). 

Note  that  one  data  card  must  exist,  so  that  when  running  the 
"standard  ;ase",  at  least  one  of  the  built-in  values  must  be 
repeated  in  the  input. 

Referring  to  Table  XXXV,  card  3  represents  data  input  for  loca¬ 
tions  5,  6,  7,  8,  9.  Since  RKML  (location  5)  is  the  only  value 
to  be  changed  from  the  built-in  data,  only  the  first  data  fiold 
is  used,  the  others  remaining  blank.  The  next  data  item  is  SBARA 
in  location  16  and  cannot  be  entered  on  the  same  card,  so  a  new 
set  of  five  items  is  inserted  on  the  next  card. 

Considerable  flexibility  exists  for  entering  input  data  on  the  same 
or  separate  data  cards.  For  instance,  RKE  started  with  the  built-in 
value  of  4.0,  was  changed  to  1.0  by  card  19,  changed  back  to 
4.0  by  card  23  and  finally  changed  back  to  1 .0  by  card  29. 
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APPENDIX  IV 

EXAMPLES  OF  USE  OF  PROGRAM 


A4. 1  Summery 

Two  groups  of  examples  are  given  in  this  Appendix,  to  illustrate  the  input 
ond  output  options  available.  The  first  group  uses  the  standard  data  and 
the  second  group  uses  realistic  C-141  gust  load  data  with  realistic  strength 
data.  One  of  the  cases  in  this  second  group  was  arranged  to  hove  com¬ 
parable  data  to  that  required  by  the  original  program  of  reference  1,  and 
the  corresponding  case  was  run  using  that  program. 

A4.2  Standard  Data 

The  input  data  is  shown  in  Table  XXXVI  .Three  cases  were  run  to  illustrate 
the  output  options.  Case  1  makes  no  changes  to  the  built-in  data  (element 
]  is  repeated  to  provide  DECRD  with  input);  the  output  tables  commence  at 
line  13,  the  highest  value  of  X  with  PXL  of  1.0,  and  ends  at  the  line  where 
PRS  reaches  1.0.  The  last  line  is  also  printed. 

Case  2  cells  for  every  line  to  be  printed  by  a  negative  value  for  RNB 
(element  4);  plotted  output  is  requested  by  1.0  at  location  140.  The 
third  case  restores  element  140  to  zero  (no  plots),  but  inputs  DX  as  nega¬ 
tive  to  call  for  the  short  output  with  the  data  and  tables  omitted. 

The  ot<>put  appears  in  Tables  XXXVII  through  XXXIX  and  the  plots  of 
case  2  are  given  in  Figure  99. 

A4.3  C-141  Examples 

The  input  data  is  that  of  Tablo  XXXV.  Eight  cases  were  run  in  sequence, 
the  first  four  having  three  tests  surviving  1 .5  times  the  unfactored  load 
and  the  last  four  having  three  test  failures  at  that  load.  The  error  furetion 
for  cases  1-3  and  5-7  was  a  double-family  distribution,  the  means,  co¬ 
efficients  of  variation  and  contributions  of  the  two  families  being  varied. 
Cases  4  and  8  employ  the  standard  Jablecki  function  of  reference  1. 

The  output  of  cases  1 ,  4  and  8  is  given  in  Tables  XL,  XU  and  XLII . 

The  plotted  outputs  of  cases  1  and  4  are  shown  in  Figure  100  and  101. 

Table  XLII  I  summarizes  the  output  of  all  eight  cases,  together  with  the 
output  from  two  runs  of  the  original  program.  These  used  a  Weibull  strength 
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distribution  (skewed  towards  tower  strength),  but  whereas  case  9  used  a 
Weibull  toad  distribution  (skewed  towards  lower  loads),  case  10  used 
a  log-normal  load  distribution  (skewed  towards  higher  loads)  as  being 
more  representative  of  the  data  used  for  Case  4. 

Examination  of  the  values  shows  that  in  spite  of  the  differences  in  data, 
in  error  functions  and  in  the  two  programs,  relatively  little  differences 
exist  in  the  reliabilities  "demonstrated"  by  the  test  results. 
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(a)  With  probable  discrepancy,  no  test 
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APPENDIX  V 

LOAD  AND  STRENGTH  DATA 


A5. 1  Introduction 

This  Appendix  gives  examples  of  the  types  of  data  available  for  the 
derivation  of  loads  and  strength  distributions,  and  the  way  in  which 
suitable  simple  statistical  equations  can  be  fitted  to  such  data  collec¬ 
tions.  The  double-family  techniques  of  Appendix  I  have  been  used 
for  the  examples,  with  the  specific  equation  of  the  Gumbel  distribution 
(the  first  asymptotic  theory  of  extremes). 

For  the  load  distributions,  the  distribution  of  maximum  values  has  been 
chosen  as  giving  better  representation  of  the  most  significant  region  -  the 
high  load  end  of  the  distribution.  In  the  case  of  the  strength  data,  the 
most  significant  area  is  that  of  low  strength  and  the  distribution  of  minimum 
values  has  been  selected.  The  logic  of  this  choice  is  obvious:  failure  is 
most  likely  to  be  the  result  of  a  high  load  or  a  low  strength;  the  exact 
representation  of  the  low  loads  and  high  strengths  will  not  contribute  to 
abetter  assessment  of  the  reliability. 

Irrespective  of  any  formal  mathematical  arguments,  it  is  essentiol  that  the 
equations  used  do  provide  a  "reasonable"  fit  to  the  data  Gumptions  of 
particular  distributions  without  convincing  evidence  of  tn.ir  validity  con 
only  lead  to  repetition  of  Disraelis  famous  criticism.* 

A5.2  Load  Data 

a.  Gust  and  maneuver  load  distributions  for  the  C-141 A  have  been 
described  in  Section  X.  Data  for  landing  impact,  taxi,  take¬ 
off  and  landing  run-out  are  described  in  this  Appendix.  All  four 
load  conditions  have  been  used  for  main  landing  gear  loads,  and 
the  three  ground  conditions  hove  been  used  to  derive  wing  root 
bending  moments. 


*w  There  are  three  kinds  of  lies:  lies,  damned  lies  and  statistics." 


b.  Landing  Impact 

Sink  rate  data  from  5345  Ground  Loads  Survey  landings  have 
been  assembled  and  fitted  by  a  double-family  Gumbel  distri¬ 
bution  as  shown  in  Figure  102. 

These  results  have  then  been  combined  with  fuel  and  cargo 
data  from  the  usage  analysis  to  determine  the  vertical  and 
drag  loads  on  the  main  landing  gear,  the  resulting  distributions 
being  presented  in  Figure  103.  The  design  lifetime  of  12000 
landings  is  used  as  the  return  period.  The  original  basic  design 
case  for  the  gear  was  a  10  ft/sec  landing  at  design  landing  weight 
(257,100  lb.).  This  limit  load  and  the  design  ultimate  (1.5  times 
limit)  load  are  shown  for  comparison  wiln  limit  and  omega  con¬ 
ditions  defined  at  the  suggested  probabilities  of  10°  and  10’^  per 
lifetime. 


c. 


It  can  be  seen  that  the  original  design  limit  and  ultimate  conditions 

3 

hove  observed  probabilities  of  approximately  10"°  and  10""  per  life¬ 
time,  so  that  the  present  design  limit  condition  approximates  the 
suggested  omega  (overload)  condition,  and  the  present  ultimate  con¬ 
dition  exceeds  the  omega  condition  by  50  per  cent. 

Taxi,  Takeoff  and  Runout 


An  arbitrary  2,0g  static  taxi  requirement  provided  the  design  down  - 
bending  case  for  the  C-141  inner  wing  and  landing  gear  vertical 
load.  Power  spectral  analyses  were  conducted  to  assess  llie  probability 
of  such  a  condition,  using  the  methods  of  reference  16.  Four  types 
of  surfacs,  were  assumed,  ranging  from  "prepared-smooth"  to  "unpre¬ 
pared-rough".  A  20  knot  taxi  speed  was  assumed;  the  takeoff  and 
runout  velocities  were  varied  with  weight,  and  the  takeoff  analysts 
Included  appropriate  lift. 


The  usage  data  were  then  employed  to  solve  for  the  wing  and  gear  load 
spectra,  using  the  following  exceedence  equations: 


N(y)  -  Nq  •  T  •  \  PRW  •  exp 
v  RW=1 
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Cumulative  Probability  of  Sink  Rote  >  A  Given  Value  Per  anding 


PIG Ufc£  \02  C-I4I  LANDING  IMPACT  SINK'  SATE  SPECTRUM 
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PftefedbfSHy  <s?  lyood  Si  A  Given  Value  In  12,000  Lendings 


Take-off  and  run-out: 

4 

N(y)  =  Y  P  •  N  *  T 


Where: 


v=! 

N(y)  = 


N 


RW 


r  4 

V 

L*u 

RW-1 


RW 


exp 


peak  spectrum  for  a  particular  data  block 
(fuel-cargo  combination) 

®  characteristic  frequency  of  aircraft  response 
for  a  given  data  block  and  velocity 

=  total  time  in  a  particular  data  block 

=  fractional  time  for  each  of  the  four  runway 
roughness  levels 

=  fractional  time  for  a  particular  data  block/ 
velocity  combination 

*  total  peak  load 


A5.2 


t  .Gg  static  load  for  a  particular  data  block 


yv  °  l.Og  mean  load  for  a  particular  data  block 
V  and  velocity 

^RW  3  r*m.s.  incremental  lead  response  variable 

far  o  given  dato  block,  velocity  and  runway 
roughness  level . 

Figures  104  and  505  shew  the  resultant  wing  downibending  ;£>cc*fum 
and  the  main  gear  vertical  loud  spectrum,  respectively,  toads  corre¬ 
sponding  to  the  suggested  limit  and  omega  levels  ore  indicated,  the 
forager  being  less  than  hvo- thirds  of  the  original  2.0  static  taxi  loads 
while  the  latter  are  less  than  three-guarters  of  these  design  loads. 


It  is  also  seen  that  the  original  design  limit  load  hot  o  probability 

7 

of  exuv'edonce  of  less  than  10"  per  12000  landing  lifetime,  sug¬ 
gesting  that  the  2,Gg  static  taxi  case  isve*y  conservative  compared 
with  probabilistic  taxi,  taxe-off  u<  J  landing  run -cut  loads. 


*f  spectral  density  analyses  are  used  to  derive  design  conditions,  care 
must  be  taken  to  ensure  shat  the  basic  parameters  do  represent  all 
passible  sources  of  (coding,  including  such  events  as  tawing  over 
damaged  surfaces,  etc. 
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%  Limit  MLG  Vertical  Load 


FIGURE  105  C— 141  MAIN  LANDING  GEAR  VERTICAL  LOAD  PROBABILITY  SPECTRA 

FOR  TAXI,  TAKEOFF  AND  RUNOUT 


A5.3  Mqteriol  Strength 


a.  Four  typical  sets  of  test  results  have  been  extracted  and  subjected 
to  analysis  to  derive  equations  for  double-family  Gumbel  distri¬ 
butions. 

b.  Aluminum  Alloy  7079-16 

Figure  106  (o)  shows  the  frequency  distribution  of  the  183 
observations,  with  the  corresponding  probabilities  of  a  lower 
value  in  Figure  106(b),  using  the  transformed  parameter,  Y, 
to  obtain  a  plot  on  the  Gumbel  paper  scale.  The  fitted  single¬ 
family  distribution  is  represented  by  the  straight  line  on  this 
latter  figure,  and  is  also  shown  on  the  frequency  distribution 
figure. 

The  double-family  distribution,  shown  on  both  figures,  was 
derived  from  the  Lockheed-Georgia  Company  program  EVDIS 
and  shows  better  correlation  in  the  significant  region  of  the 
low-strength  tail .  An  even  better  fit  could  easily  be  obtained 
by  some  further  adjustment  of  the  weaker  family. 

c.  300  VAR  Steel,  280  KSI 

Figures  107  show  similar  improvements  over  the  use  of  a  single 
family  distribution.  The  fitted  double  family  has  the  smaller 
sub-family  at  the  high  strength  end. 

d.  Titanium  Sheet  at  80°F 

Figures  108  contain  corresponding  functions  for  this  material. 

e.  Boron  composite  specimens  in  longitudinal  flexure 

A  group  of  68  test  results,  from  specimens  fabricated  over  a  period 
of  almost  one  year,  was  assembled  and  analyzed  in  order  to  compare  the 
characteristics  with  those  of  typical  metallic  materials.  Figure  109(a) 
shows  the  fitted  normal  distribution;  this  is  an  excellent  fit  in  the 
vicinity  of  the  mode,  but  misses  the  three  lowest  values  completely, 
as  shown  in  the  lower  figure.  Figure  109(b)  illustrates  the  better  fit 
obtained  with  a  (skewed)  Gumbel  distribution  of  minima.  Two  types 
of  double-family  distribution  were  then  tried;  two  families  added,  ns 
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FIGURE  107  300  VAR  STEEL  STRENGTH 
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(b)  CUMULATIVE  DISTRIBUTION 


FIGURE  107  (CONCLUDED) 


(a)  FREQUENCY  DISTRIBUTION 


FIGURE  108  TITANIUM  SHEET  STRENGTH 


shown  in  figure  109(c) ,  results  in  no  improvement;  with  family 
B  subtracted  from  Family  A,  however,  as  illustrated  by  figure 
109(d)  ,  a  better  fit  is  obtained. 

g.  Table  XUV  summarizes  the  results  of  the  study  of  material  strength 
date.  The  99  per  cent  exceedence  values  (ignoring  the  confidence 
level)  are  shown  for  comparison  purposes.  In  the  case  of  the  boron 
composite,  whose  skewness  is  the  most  pronounced,  it  can  be  seen 
that  the  assumption  of  normality  could  lead  to  a  design  value  which 
is  significantly  greater  than  that  derived  from  a  skewed  distribution. 

A5.4  Joint  Strength  Data 

a.  A  series  of  test  data  sets  was  examined  as  a  possible  approach  to  the 
selection  of  a  fabrication  variation.  In  many  instances,  groups  of 
riveted  Joint  specimens  will  be  made  from  material  of  a  single  batch, 
so  that  relatively  little  material  strength  scatter  could  be  anticipated. 

b.  AD5  rivet,  t  -  0.05  inches 

A  set  of  twer^y  test  results  on  riveted  joints  using  AD5  rivets  in  0.05 
inch  7075-T6  sheet  was  analyzed  with  the  results  shown  in  Figure  110. 

c.  06  rivet,  t  =  .09  inches 

A  second  group  of  ten  riveted  joint  results  was  also  analyzed,  figure 
1 1 1  showing  the  ability  of  the  double-family  method  to  represent 
distributions  of  a  distinctly  bi-modal  character;  the  number  of  data 
points  is  too  small  for  definitive  results. 

d.  Taper-lok  fasteners 

A  series  of  groups  of  tests  was  examined,  the  results  in  each  group 
(generally  of  ten  specimens)  being  expressed  as  a  fraction  of  the 
group  mean.  The  results,  in  figure  112  ,  clearly  show  the  unusual 
distribution  and  the  way  in  which  the  tail  Is  reproduced  by  a  double¬ 
family  distribution  with  the  second  family  subtracted  rom  the  first, 

e.  Lockbolts 

Figures  113  ar>d  114  give  the  results  of  analyses  of  two  groups 
of  lock-bolt  joint  tests;  each  containing  ten  specimens  only.  While 
inconclusive  as  definitive  values,  the  use  of  the  two  types  of  double- 
family  distribution  is  demonstrated. 
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68  SPECIMENS 
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(c)  CUMULATIVE  PROBABILITY 


FIGURE  U2  (CONCLUDED) 


10  Observations 


A5.5  Conclusions 

a.  The  inaccuracies  which  can  result  from  the  assumption  that  a  data 
sample  has  any  particular  distribution  shape  must  be  emphasized; 
selection  of  a  distribution  must  be  based  on  examination  of  the 
data  to  be  represented. 

The  importance  of  representing  the  tails  of  the  distributions  (upper 
end  for  ionds,  lower  end  for  strength)  needs  special  care  in  the 
prediction  of  failure  risks,  rince  it  is  these  toils  which  are  most 
important.  • 

b,  ;  ihe  realization,  that  the  equations  are  simply  a  means  of  expressing 
:  the  characteristics  of  the  important  parts  of  observed  distributions  in 

a  convenient  algebraic  fashion  (whether  the  equations  are  normal, 
log-normal,  V/eibulls  Gumbei,  heorson.  en-.  is ef  no  consaqusnce) 
will  avoid. the  cho?^©  made  by  Andrew  on  the  mon  who  “usat  ■ 
statistics  as fer  support,  no*  fa.-  •>  ’ 
illumination*" 
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